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Abstract: We analyze the production of Higgs particles at the early stage of the CERN 

■ large Hadron Collider with a 7 TeV center of mass energy (IHC). We first consider the 
case of the Standard Model Higgs boson that is mainly produced in the gluon-gluon fusion 
channel and to be detected in its decays into electroweak gauge bosons, gg ^ H ^ 
WW, ZZ, 77. The production cross sections at y/s = 7 TeV and the decay branching 
ratios, including all relevant higher order QCD and electroweak corrections, are evaluated. 
An emphasis is put on the various theoretical uncertainties that affect the production 
rates: the significant uncertainties from scale variation and from the parametrization of 
the parton distribution functions as well as the uncertainties which arise due to the use of 
an effective field theory in the calculation of the next~to-next-to-leading order corrections. 
The parametric uncertainties stemming from the values of the strong coupling constant and 
the heavy quark masses in the Higgs decay branching ratios, which turn out to be non- 
negligible, are also discussed. The implications for different center of mass energies of the 
proton collider, ^/s = 8-10 TeV as well as for the design energy ^/s = 14 TeV, are briefly 
summarized. We then discuss the production of the neutral Higgs particles of the Minimal 
Supersymmetric extension of the Standard Model in the two main channels: gluon-gluon 
and bottom quark fusion leading to Higgs bosons which subsequently decay into tau lepton 
or 6-quark pairs, gg,bb — )• Higgs — )• T~^T~,hh. The Higgs production cross sections at the 
IHC and the decay branching ratios are analyzed. The associated theoretical uncertainties 
are found to be rather large and will have a significant impact on the parameter space of 
the model that can be probed. 
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1. Introduction 

The Standard Model (SM) of the electroweak and strong interactions crucially relies on the 
Higgs mechanism in order to spontaneously break the electroweak symmetry and generate 
in a gauge invariant way the elementary particle masses [1, 2]. In this minimal realization 
of electroweak symmetry breaking, one complex Higgs doublet field is introduced which 
implies the existence of a single neutral scalar particle, the Higgs boson H. In supersym- 
metric theories [3], that are widely considered as the most attractive extensions of the 
SM as they protect the Higgs boson mass against large radiative corrections and stabilize 
the hierarchy between the electroweak and Planck scales, the Higgs sector is extended to 
contain at least two Higgs doublet fields. The Minimal Supersymmetric Standard Model 
(MSSM) predicts the existence of five Higgs particles: two CP-even Higgs particles h and 
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H, a CP-odd A boson and two charged H particles [4-7]. The search for these new scalar 
particles is the main goal of present high-energy colliders. 

With its very successful operation in the last years, the Fermilab Tevatron pp collider 
has now collected a substantial amount of integrated luminosity which allows the CDF 
and DO experiments to be sensitive to Higgs particles. Exclusion limits beyond the well 
established LEP bounds [8] are now being set on Higgs masses both in the SM [9] and in 
the MSSM [10]. The CERN proton-proton collider successfully started operation but at 
a reduced center of mass energy of 7 TeV [11]. Not to confuse it with the LHC, which 
is expected to operate at the design energy close to 14 TeV, we will call this early stage 
machine the IHC for large (or littler) Hadron Collider. The IHC will also be sensitive 
to the SM Higgs particle and will start to be competitive with the Tevatron once it will 
accumulate the expected 1 fb^"*^ of data [11-14]. However, with the present expectations, 
the Tevatron and presumably also the IHC will only be able to exclude the existence of 
the SM Higgs particle in some given mass range. An undebatable discovery in the entire 
Higgs mass range favored by theoretical considerations, 115 GeV < Mh ^ 0{1 TeV), has 
to await for the full-fledged LHC with a center of mass energy of 14 TeV and 0(10 fb~^) of 
integrated luminosity. Observing the MSSM Higgs bosons at the Tevatron and the IHC will 
be more likely as, in some areas of the supersymmetric parameter space, the production 
cross sections arc much higher than in the SM [5]. 

The exclusion of Higgs mass regions relies crucially on the theoretical predictions for the 
production cross sections for the Higgs signal as well as for the relevant SM backgrounds. At 
the Tevatron, the two main search channels for the SM Higgs boson are the top and bottom 
quark loops mediated gluon-gluon fusion mechanism gg ^ H [15] with the Higgs decaying 
into WW pairs which lead to Iflu (with £ = e, ji) final states [16] and the Higgs-strahlung 
processes qq ^ VH (with V = W, Z) [17] with the subsequent H ^ bb and V ^ £ + X 
decays of the Higgs and the associated gauge bosons. Analyses performed by the CDF 
and DO experiments, under some assumptions for the production cross sections and their 
associated uncertainties, have recently excluded the Higgs mass range Mh = 158-175 GeV 
at the 95% confidence level (CL) [9]. In this mass range, the Higgs signal is mainly due to 
the gg^H -^WW -^ivlu production and decay channels. 

At the IHC, the Higgs-strahlung processes, as well as as other production channels such 
as weak vector-boson fusion and associated Higgs production with top quark pairs, have too 
small cross sections and/or are plagued with too large QCD backgrounds. Thus, in practice, 
only the gluon-gluon fusion process with the Higgs boson decaying into H^WW^Uvv, 
H ^ ZZ ^2i+X, where X stands for charged leptons, neutrinos and eventually also jets 
including 6-quark jets, and to a lesser extent — )■ 77 final states will be mostly relevant. 
At a center of mass energy ^/s = 7 TeV and with 1 fb~^ of data, recent studies by the 
ATLAS and CMS collaborations have shown that the mass range Mh fa 150-190 GeV can 
be excluded at 95% CL if no Higgs signal is observed [12, 13]. 

Nevertheless, it is well known that the production cross sections at hadron colliders 
as well as the associated kinematical distributions are generally affected by various and 
possibly large uncertainties. In a recent study of SM Higgs production at the Tevatron 
[18], it has been re-emphasized that while the signal production cross sections are well 
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under control in the Higgs-strahlung processes, the theoretical uncertainties are rather 
large in the case of the gluon-gluon fusion channel. These uncertainties are stemming 
mainly from the variation of the energy scale at which the process is evaluated (a variation 
that provides a hint of the not yet calculated perturbative higher order corrections) and 
the parametrization of the parton distribution functions or PDFs (in particular, the gluon 
density at moderate to high Bjorken-x values). At the Tevatron, these uncertainties can 
shift the central values of the gg ^ H cross section [19-22], by more than « 20% in both 
the case of the scale and the PDF uncertainty. Arguing that these are two theoretical 
uncertainties that have in principle no statistical ground and should thus not be added in 
quadrature, we have proposed a more adequate procedure to combine them which led to 
a total theoretical uncertainty of about ±40% on the gg ^ H ^ £i/£i/ signal. This overall 
uncertainty is much larger than the uncertainty assumed in the CDF and DO combined 
analysis that excluded the SM Higgs boson at the 95% CL in the mass range Mjj = 158-175 
GeV [9]. To our opinion, this exclusion limit should thus be reconsidered. 

In the present paper we extend the analysis of Ref. [18] to the case of the IHC, i.e. for 
proton-proton collisions with a center of mass energy of 7 TeV. We first update the Higgs 
production cross sections in the gg ^ H channel, including all the relevant higher order 
QCD and electroweak corrections and then analyze the various uncertainties that affect 
them. We show that the scale and PDF uncertainties, as well as the non-ncgligiblc uncer- 
tainty due to the use of an effective approach in the calculation of QCD and electroweak 
higher order corrections beyond next-to-leading order, add up to 25-30% depending on 
the considered Higgs mass range. The total uncertainty is significantly smaller than that 
obtained at the Tevatron, as a result of smaller QCD radiative corrections and a better 
knowledge of the gluon distribution function at the energies relevant at the IHC. 

We then discuss an additional source of theoretical uncertainties that has not been 
considered neither in the experimental analyses [9, 12, 13] nor in Ref. [18]: the one affecting 
the Higgs decay branching ratios. Indeed, while the Higgs decays into lepton and gauge 
boson pairs are well under control (as mainly small electroweak effects are involved), the 
partial decays widths into quark pairs and ghions arc plagued with uncertainties that are 
mainly due to the imperfect knowledge of the bottom and charm quark masses and the 
value of the strong coupling constant Ug- Updating a previous analysis [23], we show that 
at least in the intermediate mass range, Mh ~ 120-150 GeV, where the SM Higgs decay 
rates into hh and W^W~ final states have the same order of magnitude, the parametric 
uncertainties on these two main Higgs decay branching ratios are non-negligible, being of 
the order of 3 to 10% at the la level. 

For completeness, we will extend our analysis and explore the implications for SM 
Higgs searches at center of mass energies beyond ^/s = 7 TeV. Not only the full-fledged 
LHC with -y/s = 14 TeV will be considered, but also intermediate energies between -y/i = 7 
and 10 TeV that are being currently considered for a very near future, with integrated 
luminosities significantly larger than 1 fb~^ [11]. The main result that we obtain is that 
while the production cross sections for the SM Higgs particle are higher at energies ^/s = 8- 
14 TeV compared to ^/s = 7 TeV, the overall theoretical uncertainties that affect them are 
approximately the same. 
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In the MSSM, the chances of observing at the IHC the neutral Higgs particles, that 
we will collectively denote by are significantly higher than in the case of the SM Higgs 
boson. Indeed, for large values of the ratio of the vacuum expectation values of the two 
MSSM Higgs fields, tan/3 = V2/v\, the couplings of some of the Higgs bosons to isospin 
down-type particles are strongly enhanced. For not too heavy Higgs particles, this leads 
to production cross sections that are possibly orders of magnitude larger than in the SM, 
in processes in which the MSSM Higgs particles couple to bottom-quarks [5,6]. This is 
the case of the gluon-gluon fusion mechanism for neutral Higgs production gg ^ ^ which, 
in this context, dominantly proceeds through 6-quark triangular loops [15, 20], as well as 
bottom-quark fusion, in which the bottom quarks are directly taken from the protons in a 
five active flavor scheme, 65 — )■ $ [24-26]. For the relevant high values of tan /?, the $ bosons 
decay almost exclusively into hb and r+r^ final states with, respectively, ~ 90% and ~ 10% 
branching ratios. This leads to the dominant production channels, gg, bb ^ ^ t^t^ , bb. 
Again, besides the uncertainties due to scale variation and the parametrization of the PDFs, 
additional uncertainties originate from errors in the extraction of the bottom quark mass 
mb and the strong coupling constant ag. These uncertainties affect both the production 
cross sections, as the amplitudes involve the Higgs-66 Yukawa coupling oc irib, and the 
$ T^T~ , bb decay branching fractions. 

A second aim of this paper is to analyze in detail these two main MSSM Higgs pro- 
duction channels at the IHC: the production cross sections, including the relevant higher 
order radiative corrections are updated, the important QCD effects are summarized and 
the various theoretical uncertainties are evaluated. We show that the uncertainties from 
scale variation and the parametrization of the PDFs in the cross sections are at least as 
large as in the SM gg —> H case. Additional parametric uncertainties in both the produc- 
tion and the decay rates, as well as a freedom in the choice of the renormalization scheme 
that defines the 6-quark mass, lead to a total theoretical uncertainty that is even larger, 
of the order of 40 to 50% depending on the Higgs masses. This uncertainty will thus have 
a significant impact on the MSSM parameter space that can be probed at the IHC. 

The rest of the paper is organized as follows. In the next section we discuss the pro- 
duction of the SM Higgs boson at the IHC with y/s = 7 TeV, mainly concentrating on 
the dominant gg ^ H production channel; the total inclusive cross section including the 
relevant higher order contributions is updated, the various theoretical uncertainties eval- 
uated and combined, and the implications for higher center of mass energies summarized. 
In section 3, the parametric uncertainties in the SM Higgs decay branching ratios axe sum- 
marized and combined with the uncertainties in the production cross sections. In section 
4, we discuss the production of the MSSM neutral Higgs bosons: the cross sections in the 
two production channels, gg ^ ^ and 66 — )• $ are analyzed and the associated uncertain- 
ties from scale and scheme dependence, the PDF uncertainties as well as the parametric 
uncertainties are discussed; the errors in the Higgs branching ratios in the main decay 
modes $ t^t~ and ^ ^ bb are addressed and the way to combine all these theoretical 
uncertainties in the production and in the decays is discussed. Here, we only discuss the 
case of the IHC with = 7 TeV and simply summarize the implications at slightly higher 
energies, ^/s = 8-10 TeV. Finally, a conclusion is given in the last section. 
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2. The SM Higgs boson at the IHC 



2.1 The Higgs production cross sections 

The main production channel for the Higgs particles at hadron colliders is the gluon-gluon 
fusion channel, gg — > H, which proceeds through triangular heavy quark loops [15]. In 

the SM, it is dominantly mediated by the top quark loop contribution, with a bottom 
quark contribution that does not exceed the 10% level at leading order. This process is 
known to be subject to extremely large QCD radiative corrections [19-22, 27-29] that can 
be described by an associated i^T-factor defined as the ratio of the higher order (HO) to 
the lowest order (LO) cross sections, consistently evaluated with the value of the strong 
coupling as and the parton distribution functions (PDF) taken at the considered order, 

KuO = <7"°|(aHO,PDFHO) / f^^°|(aLO,pDFLO) (2.1) 

The next-to-leading-order (NLO) corrections in QCD are known both for infinite [19] and 
finite [20] loop quark masses and, at ^/s = 7 TeV, lead to a if-factor -FCnlo ~ 1.8 in the low 
Higgs mass range, if the central scale of the cross section is chosen to be the Higgs mass. 
It has been shown in Ref. [20] that working in an effective field theory (EFT) approach in 
which the top quark mass is assumed to be infinite is a very good approximation for Higgs 
mass values below the tt threshold Mh < 2mt, provided that the leading order cross section 
contains the full rrit and irih dependence. The calculation of the next-to-next-to-leading- 
order (NNLO) contribution has been done [21] only in the EFT approach Mh <C 2mt and, 
at ^/s = 7 TeV, it leads to a ~ 25% increase of the cross section-'^, i^NNLO ~ 2.5. The 
resummation of soft gluons is known up to next-to-next-to-leading-logarithm (NNLL) and, 
again, increases the cross section by slightly less than 10% [22, 28]. The effects of soft-gluon 
resumation at NNLL can be accounted for in a^^^^{gg H) by lowering the central value 
of the renormalization and factorization scales^, from /xq = Mh to /xq = \Mh- We will 
thus choose for definiteness the value 

fio = m = fiF = ^Mh (2.2) 

for the central scale of the process which, in passing, also improves the convergence of the 
perturbative series and is more appropriate to describe the kinematics of the process [33]. 
The electroweak corrections are known exactly up to NLO [34, 35] and contribute at the 

^Thc QCD corrections to gg ^ H at ^/s — 7 TeV arc thus smaller than the corresponding ones at 
the Tevatron as the /sT-factors in this case are -K^nlo ~ 2 and /Cnnlo ~ 3 (with a central scale equal to 
Mh)- At the LHC with ^fs = 14 TeV, the it'-factors are even smaller, /Tnlo <^ 1-7 and Knni,o 2. The 
perturbative series shows thus a better (converging) behavior at LHC than at Tevatron energies. 

^We will include the soft-gluon resummation contributions in this indirect way since, as discussed in 
Ref. [18], we would like to stick to the fixed order NNLO calculation for two main reasons: i) there are not 
yet parton distribution functions which include soft-gluon resummation and it appears inconsistent to fold 
a partonic cross section with PDFs that are not at the same order of perturbation theory (although the 
effects might be small in practice [30]) and n) the soft-gluon resummation is not available for the Higgs-l-jet 
production cross sections and/or the Higgs cross sections including kinematical cuts (both are known only 
at NNLO [31,32]) which, ultimately, are the basic experimental inputs. 
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Figure 1: The total cross sections for Higgs production at the IHC with ^/s = 7 TeV as a function 
of the Higgs mass. The MSTW set of PDFs has been used and the higher order corrections are 
included as discussed in the text. 

level of a few percent; there are also small mixed NNLO QCD-electroweak effects which 
have been calculated [27] in an effective approach valid for Mh ^ M^. 

Our calculation^ of the gg ^ H production cross section at the IHC, including these 
higher order corrections strictly follows the one performed in Ref. [18] for the Tevatron: the 
starting point is the Fortran code HIGLU [37] which evaluates the gg ^ H cross section at 
exact NLO in QCD (i.e. with the exact contributions of the top and bottom quark loops) to 
which we add the NNLO QCD contribution in the infinite top-quark mass limit, but with 
the LO normalisation containing the exact rrit and mi, dependence; grids for the exact NLO 
electroweak and the mixed QCD-electroweak corrections are then implemented within the 
code (they have been added in the partial factorization scheme [34], see section 2.2.2). The 
production cross sections are shown at the IHC with = 7 TeV in Fig. |l] for an updated 
value of the top quark mass and when the partonic cross section is folded with the latest 
NNLO MSTW2008 pubHc set of PDFs [38]. The renormalization and factorization scales 
are fixed to the central values = /^i?, = ^Mh- Our results agree with those given in 
Refs. [27, 28] and updated in Ref. [14] within a few percent. 

For completeness, we also display in Fig. |l]the cross sections for the three other Higgs 
production channels at hadron colliders that we evaluate using the programs of Ref. [37]: 
i) the Higgs-strahlung processes qq — )• HV with V = W,Z that are known exactly up to 

NNLO in QCD [6, 39, 40] and up to NLO for the electroweak corrections [41]; they are 

evaluated at no = Mhv (the invariant mass of the HV system) for the central scale'^; 

■^Other recent updates of the gg H cross section can be found in Ref. [27, 28, 36]. 

*Here, the QCD Jf-factors are moderate, A'nnlo ~ 1-5 and the electroweak corrections reduce the cross 
section by an amount of « 3 — 8%. For the evaluation of the cross section, we have used the NLO code V2HV 
[37] in which we implemented these higher order contributions. Our results agree with those of Ref. [14]. 
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ii) the weak vector boson fusion channel qq — t- Hqq evaluated at the scale //q = Qv (the 
momentum transfer at the gauge boson leg) in which only the NLO QCD corrections 
[6, 42] have been included; the NNLO corrections have been found to be very small 
[43] and we omit the electroweak corrections [44]; 
Hi) associated Higgs production with top quark pairs in which only the leading order order 
contribution is implemented but at a central scale /io = \ {^H + 2rn-t), which is a good 
approximation at these energies as the NLO K-ia.ctoT is very close to unity [45,46]. 

As can be seen, the gg ^ H process is by far dominating in the entire Higgs mass 
range, with a cross section that is one to two orders of magnitudes larger than in the other 
production channels. For low Higgs masses, Mh ^ 2M\\r, it leads to more than 10.000 
events for a luminosity of order 1 fb^^. Table |l| displays the values of the cross sections for 
the five Higgs production processes for a selection of Higgs masses relevant at the IHC. 

2.2 The theoretical uncertainties 

2.2.1 Higher order contributions and scale variation 

In the calculation of production cross sections and kinematical distributions at hadron col- 
liders, as the perturbative series are truncated and the results are available only at a given 
perturbative order, there is a residual dependence of the observables on the renormaliza- 
tion scale which defines the strong coupling constant Ug and on the factorization scale 
/Xi? at which the matching between the perturbative matrix elements calculation and the 
non-perturbative parton distribution functions is performed. The uncertainty due to the 
variation of these two scales is viewed as an estimate of the unknown (not yet calculated) 
higher-order terms and is rather often the dominant source of theoretical uncertainties. 

In general, starting with the median scale = M-R = /^o for which the central value 
of the cross section is obtained, the two scales fin and fip are varied within the interval 
/^o/'^ ^ ^ KXfMo with the value of the constant /t = 2,3,4, etc... to be chosen. 

The additional restriction 1/k < fip/ fiji < k is also often applied^ . The choice of the scale 
variation domain and hence the constant factor k, is rather subjective and a possibility 
would be to adopt the k value which allows the uncertainty band due to the scale variation 
of the lower order cross section to incorporate the central value of the cross section (i.e. 
with fiR = = ^q) calculated at the highest available order. 

In our analysis, we have adopted the central scale value fJ-R = fip = fJ-o = \,^H as 
discussed earlier. At ^/s = 7 TeV, for the scale uncertainty band of (y^^^{gg — )■ H) to catch 
the central value of (y^^^^{gg — t- H) for /ig = ^Mh, a value k = 2 is sufficient. Adopting 
the range ^/^o < fJ-R, fJ-F < with //q = for the cross section a^^^^{gg — )• H), 

one obtains the scale variation displayed in Fig. |2| as a function of Mh', it is compared 
to a variation with a factor k = 3. In the insert, shown are the maximal and minimal 
variations of (T^^^^{gg — s- H) compared to the central value. One sees that for k = 2, a 
scale uncertainty of ~ ±10% is obtained in the low mass range, Mh ~ 120 GeV, which 

^In the case of the gg ^ H process, the maximal (minimal) cross sections at a given fixed order in 
perturbation theory (in particular at NNLO) is obtained for the scale choices ^r, = ^lf = yLo/hi (k/iq) and, 
hence, this restriction is irrelevant in practice. 
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Table 1: The total Higgs production cross sections (in fb) in the processes gg — >■ H, vector-boson 
fusion qq — >• Hqq, Higgs-strahlung qq — >• HW.HZ and associated production pp — >• ttH at the IHC 
with y/s = 7 TeV for given Higgs mass values (in GeV) with the corresponding central scales 
described in the main text. The MSTW sets of PDFs have been used at the relevant order. 

decreases to the level of « —8%, +4% at high masses, Mh ~ 500 GeV. Note that if the 
domain for scale variation were extended to k = 3, the uncertainty would have increased 
to « ±17% in the low Mh range^ as also shown in the figure. 

®In the analysis of Ref. [18] for the case of the Tevatron, the value k = 3 was adopted instead of k = 2 
and, hence, a larger domain for the scale variation was assumed as the higher order QCD corrections in 
""(ffff H) are larger at the Tevatron compared to the IHC. 
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Figure 2: The scale uncertainty band of (T^^^'^{gg — > H) at the IHC as a function of Mh for a 
scale variation in the domain AIh /{2k) < ^r, < kx ^Mh for k ~ 2 and 3; in the insert, shown 
are the relative deviations from the central value with a scale fin. = I-I-f = \Mh- 

2.2.2 The use of an effective theory approach 

Another source of theoretical uncertainties is specific to the gluon-gluon fusion mechanism 
in which the cross section beyond the NLO approximation is evaluated in an effective 
theory (EFT) approach where the particles that are running in the loops are assumed to 
be much heavier than the produced Higgs boson. At NLO, the approximation mt ^ Mh 
for the contribution of the top quark in the loop is rather good for Higgs masses below the 
tt threshold, Mh ^ 340 GeV, in particular when the full quark mass dependence of the 
leading order cross section cf^^^^ct is taken into account [20]. At NNLO, this approximation 
for the top quark contribution seems also to be accurate as studies of the effect of a finite 
mt value in expansions of Mh / {"i^mt) have shown a difference below the percent level with 
respect to the EFT calculation for Mh < 300 GeV [47]. 

Nevertheless, the EFT approach should definitely not be valid for Higgs masses beyond 
the tt threshold, where the gg ^ H amplitude develops imaginary parts. This can be seen 
at NLO where both the exact and the approximate results are known. We will thus include 
an uncertainty related to the use of the EFT approach for Mh ^ 2mt which is taken as 
the difference between cr^aS S'^d (t^S^oo when the exact top quark mass dependence is 
included in the LO cross section and when one rescales this difference with the relative 
magnitudes of the NLO and NNLO i^T-factors, i.e. K^^S.oo/K^f}iZ- 

In the left-hand side of Fig. |3|, shown is this difference in percent for ^/s = 7 TeV. It is 
small for Higgs masses below Mh = 400 GeV, at most 2%, it increases very rapidly with 
increasing Mh values and, for Mh ^ 500 GeV, it is larger than ~ 5%. This uncertainty 
is thus not negligible and it should be accounted for (at the Tevatron, this uncertainty has 
not been discussed [18] as only Higgs masses below 200 GeV were considered). 
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Mh [GeV] Mh [GeV] 

Figure 3: The estimated uncertainties (in %) due to the use of the EFT approach in the evaluation 
of (7{gg — > H) at NNLO for the top quark loop contribution for Higgs masses beyond 2mt (left) 
and for the bottom-quark loop contributions (right). 

Furthermore, the EFT approach is certainly also not valid for the 6-quark loop and, for 
instance, the omission of this contribution (at LO) leads to a ~ 10% difference compared to 
the exact case. Furthermore, the NLO iC-factor for the 6-quark loop, -f^b'-^io^ ~ 1.2-1.4, is 
much smaller than that for the top-quark loop. In fact, the problem arises mainly because 
of the significant negative interference between the top and bottom loop contributions (the 
bottom quark contribution itself is rather small) up to NLO and the fact that the bottom 
contribution has been factorized out in the LO cross section. Thus, including the top quark 
contribution only using the EFT approach might overestimate the NNLO correction as this 
interference component is missing (and, since it is resulting from an approximation it is 
not, in principle, accounted for by e.g. the scale variation that is performed at NNLO). 

In order to estimate the uncertainty due the missing 6-loop contribution at NNLO, we 
simply follow the previous procedure for the top quark: we assign an error on the NNLO 
QCD result which is approximately the difference between the exact result cr^act s-'^d the 
approximate result cr^^^ obtained at NLO but rescaled with the relative magnitude of 
the ET-factors that one obtains at NLO and NNLO for the top-loop, i.e. K^^^/K^}^%. 
This leads to the uncertainty on the cr^'^^^{gg — t- H) at ^/s = 7 TeV that is shown in the 
right-hand side of Fig. ^ It ranges from ~ ±3% for low Higgs values Mu ~ 120 GeV where 
the 6-loop is relatively important to less than ~ ±1% for Higgs masses above Mh ~ 300 
GeV. Note that this correction is larger than that obtained at the Tevatron [18] where a 
different central scale, = Mh-, has been adopted for its evaluation. 

In addition, there is some freedom in the choice of the renormalization scheme for the 
6-quark mass in the gg ^ H amplitude: for instance, one can use the on-shell scheme in 
which the pole mass is m^y ~ 4.7 GeV or the MS scheme in which the mass frihijnh) ~ 4.2 
GeV is adopted. This leads to a difference of ^ 1% in the 6-quark loop contribution at 
NLO that we will take as an additional uncertainty due to the scheme dependence (this 
scheme dependence will be discussed in more details in the MSSM case). 
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Finally, in the corrections to the gg ^ H amplitude, the mixed QCD-electroweak 
corrections have been calculated at NNLO [27] in an EFT approach with My/jz ^ Mh, a 
limit that is obviously not valid in practice. We thus should be cautious when using this 
correction and assign an uncertainty which is of the same size as the contribution itself [18] . 
This uncertainty is comparable in size to the difference between the electroweak correction 
calculated exactly at NLO [34] evaluated in the partial factorization scheme (where the 
correction (t^^Aew is simply added to the QCD corrected cross section) compared to the 
result in the complete factorization scheme (where the NNLO cross section is multiplied by 
1 + Aew)- This generates an additional uncertainty of f« 3% at most, which is not shown 
as it is the same as the one discussed in Ref. [18] for the Tevatron. 




2 I ^ ^ ^ ^ ^ ^ ^ 1 

150 200 250 300 350 400 450 500 
Mh [GeV] 



Figure 4: The estimated total uncertainty (in %) at IHC energies and as a function of Mh from 
the use of the EFT approach for the calculation of the gg —?' H amplitude at NNLO and from the 
scheme dependence of the 6-quark mass. 

Adding the uncertainties from these three sources linearly, the resulting overall uncer- 
tainty is displayed in Fig. Q as a function of Mjj- It amounts to ~ 7% in the low Higgs 
mass range, drops to the level of ~ 4% in the mass range Mh ~ 200-400 GeV, and then 
increases to reach the level of 6% at 600 GeV, as a result of the bad Mh ^ 2mt ap- 
proximation for the top loop contribution. Since this is a pure theoretical uncertainty, it 
should be added linearly to the uncertainty from scale variation. 

2.2.3 The PDFs and uncertainties 

One of the most important sources of theoretical uncertainties in the gg ^ H production 
cross section is the one stemming from the parametrization of the gluon densities. There 
are two quite different ways to estimate this uncertainty. 

The first one is the so-called Hessian method where, besides the best fit PDF with 
which the central values of the cross sections are evaluated, a set of 2A^pdf PDF param- 
eterizations is provided that reflect the it la variation of all (Np) parameters that enter 
into the global fit. These uncertainties are thus mostly due to the experimental errors in 
the various data that are used in the fits. 
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If one takes the latest NNLO public set provided by the MSTW collaboration [38, 48] 
for instance, the PDF uncertainty at the 90% CL that is obtained at the IHC for the NNLO 
gg ^ H cross section is shown (by the red lines and red band) in Fig. |5| as a function of 
Mh- For low Higgs masses, the uncertainty is at the level of 5% but it increases to reach 
the level of ~ 8% at high Higgs masses, Mh ~ 500 GeV. If only the (more optimistic) 
68%CL errors are to be considered, the previous numbers have to be divided by w 1.6. 




150 200 300 400 500 

Mh [GeV] 



Figure 5: The central values and the 90% CL PDF, POF+A'^'^Pa^ and POF+A'^'^Pa^ + A^^a^ 
uncertainty bands in (J^^^'~'{gg H + X) at the IHC when evaluated within the MSTW scheme. In 
the inserts, shown are the same but with the cross sections normalized to the central cross section. 

The previous Hessian errors do not account for the theoretical assumptions that enter 
into the parametrization of the PDFs and which may come from many different sources. 
These assumptions lead to an uncertainty in the parametrization that is difficult to quantify 
in a given scheme, say the MSTW scheme. A way to access this theoretical uncertainty is to 
compare the results for the central values of the cross section (and hence, using the best fit 
PDFs) when using different sets of PDFs which involve, in principle, different assumptions. 
In Fig. ^, we display the values of (7^^^'~'{gg — H) that are obtained when folding the 
partonic cross section with the gluon densities that are predicted by the four PDF sets^ 
that have parameterizations at NNLO: MSTW [38], JR [51], ABKM [52] and HERAPDF 
[53]. In the later case, two sets are provided: one with the value = 0.1176 that is close 
to the world average value [54] and the MSTW best-fit as = 0.1171 [38], and one with the 
value as = 0.1145 used also in the ABKM set [52], that is close to the preferred values that 
one obtains using deep-inelastic scattering data alone. 

As can be seen, while the differences in the cross sections are moderate in the low 
Higgs mass range, being of the order of 10% or less, there is a significant difference at 

''As we make the choice of using parton cross sections and PDFs that are consistently taken at the same 
perturbative order, we will not consider two other parameterizations, CTEQ [49] and NNPDF [50], which 
provide PDF sets only at NLO. 
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higher Higgs masses, 0(25%) for Mh ~ 500 GeV. Even with the same input or best-fit 
values, HERAPDF with = 0.1176 and MSTW results are strikingly different. For the 
large Mh values, this discrepancy is mainly due to the gluon densities at moderate to high 
Bjorken-x values which are less constrained by the data (in particular if the Tevatron high 
Et jet data are not included as is the case of the ABKM and HERAPDF sets). 
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Figure 6: The central values of the NNLO cross section cri^gg H + X) at the IHC as a function 
of Mh when evaluated in the four schemes which provide NNLO PDFs. In the inserts, shown are 
the same but with the cross sections normalized to the MSTW central cross section. 

We thus conclude when comparing Fig. ^ and Fig. |^ that, depending on how the 
uncertainties are estimated, through the 2x20 grids given in the MSTW scheme (other PDF 
sets such as ABKM for instance give approximately the same error) or by comparing the 
cross sections evaluated using PDFs from different schemes, a rather different uncertainty 
is obtained. To partly reconcile these two different results, we will follow the analysis made 
in Ref. [18] and include the additional errors on the coupling a^. 

As a matter of fact, in addition to the PDF uncertainties, one should also consider the 
ones coming from the uncertainties in the determination of the value of the strong coupling 
constant^ which, in general, is fitted together with the PDFs and this point is now well 
taken by the experimental collaborations. Indeed, since already at leading order one has 
^^^^{99 ~^ H) oc 0{a1) and the i^T-factors are large, the uncertainty in the determination 
of as may induce a non-negligible error on the final cross section. In the MSTW scheme, 
the value of as together with its experimental accuracy, is [48] 

a,(M|) = 0.1171 toZlt (68% CL) or (90% CL) at NNLO (2.3) 

*One can see that the input value of Qs is very important by comparing the two HERAPDF results in 
Fig. ^. Nevertheless, the difference in as is not sufficient to explain the behavior of the cross sections in 
the MSTW and HERAPDF cases for instance, and there is some part that is presumably also due to the 
different shapes of the gluon densities. 
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We have evaluated the 90% CL correlated PDF+A'^^Pas uncertainties using the new set-up 
provided by the MSTW collaboration [48] and the result is shown in Fig. ^ (green band 
and green lines) as a function of Mh- As one can see, the obtained uncertainties are much 
larger than the pure PDF uncertainties alone. This is particularly true at low Higgs masses 
where the uncertainty is doubled. Nevertheless, even with this additional uncertainty, one 
still cannot reconcile the MSTW and the ABKM/HERAPDF predictions for instance, in 
particular at high Higgs masses. 

However, it has been argued in Ref. [18] that one also needs to consider the theoretical 
uncertainty on that is due to the truncation of the perturbative series, the different 
heavy flavor schemes, etc... This is particularly true as one has a significant difference 
between the ag world average value and the value obtained from deep-inelastic scattering 
data alone, a difference which could be reduced if this theory uncertainty is included. The 
MSTW collaboration estimates the theoretical uncertainty to be A^^a^ = 0.003 at NLO 
[48], which gives A^^as = 0.002 at most at NNLO. Using a fixed as NNLO grid with 
central PDFs given by the MSTW collaboration, with as values different from the best- 
fit value (in the range 0.107-0.127 with steps of 0.001 and which thus include the values 
UsiM"^) = 0.1171 ± 0.002 at NNLO), we have evaluated this uncertainty. Adding it in 
quadrature to the PDF+A^^Pa^ uncertainty, we obtain a total POF+A'^^Pa^ + A*^a<j 
uncertainty of ~ 11% — 15% depending on the Higgs mass. At least for not too heavy Higgs 
bosons, this larger uncertainty reconciles the MSTW and ABKM/HERAPDF predictions^. 

Note that for Mjj < 200 GeV, we obtain smaller uncertainties compared to the ~ ±15- 
20% uncertainty of o'^^^'^{gg — t- H) at the Tevatron, a reflection of the better control on 
the behavior of the gluon density at moderate-x values (relevant for the IHC) compared 
to high-x values (relevant for the Tevatron). 

We note also that the result that we obtain here for the RDF+A^^Pa^ + A^'^a^ uncer- 
tainty is larger than what one obtains using the recommendation of the PDF4LHC group 
[56], that is, to take the envelope of the PDF+A°^Pas values provided by the MSTW, 
CTEQ and NNPDF sets or more simply (and more consistently at NNLO), take the 68%CL 
MSTW PDE+A^^^Pa^ error and multiply it by a factor of two. To our opinion, this estimate 
of the uncertainty is rather optimistic and, as can be seen from Fig. ^, cannot account for 
the difference between the predictions using the various NNLO parameterizations. 

In the rest of our analysis, we will take the 90% CL PDF+A'^^Pas + A^'^a^ uncertainty 
evaluated as above as the total PDF uncertainty in the determination of the gg ^ H 
production cross section. However, we will view this uncertainty as a measure of the 
difference between the various possible parameterizations of the PDFs and, hence, consider 
it as a true theoretical uncertainty with no statistical ground as will be discussed shortly. 

®Note that we could also take into account the effects due to the bottom and charm quark masses on 
the PDFs. Indeed, a change in the fitted masses for these quarks may affect the gluon splitting and hence 
alter the shape of the gluon-gluon luminosity in gg — >■ H. We have estimated quantitatively this effect by 
using the mt and dependent MSTW PDFs [55] with = 1.40 ± 0.15 GeV and = 4.75 ± 0.05 GeV. 
In the case of the charm quark, we obtain a ~ 0.2% change at Mh ~ 115 GeV and at most « 1.5% change 
at Mh ~ 500 GeV. In the case of the 6-quark, the change is below the percent level (see also section 4.3). 
These uncertainties can be neglected as we adopted the attitude of adding them in quadrature with the 
ones due to the PDF and Us which are much larger. 
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2.3 The combined uncertainty 

A very important issue that remains is the way to combine the various theoretical uncer- 
tainties on the cross section discussed in the previous subsection. Let us first reiterate an 
important comment: the uncertainties associated to the PDFs in a given scheme should be 
viewed as purely theoretical uncertainties despite of the fact that they are presented as the 
la or more departure from the central values of the data included in the PDF fits. Indeed, 
this uncertainty should be interpreted as being due to the parametrization of the quark and 
gluon densities and, hence, to the theoretical assumptions in determining these densities 
and equivalent to the spread that one observes when comparing different parameteriza- 
tions. Thus, the PDF uncertainties should be considered as having no statistical ground^*^ 
and they cannot be combined in quadrature with the two other purely theoretical errors, 
namely the scale and the scheme (EFT) uncertainties. This should be contrasted to what 
has been done, for instance, by the CDF and DO collaborations in their combined Higgs 
analysis at the Tevatron [9], where the scale and PDF uncertainties have been added in 
quadrature. 

Here, we discuss three possible ways to combine the various uncertainties. 

A) A reasonable procedure to combine the uncertainty from scale and scheme dependence 
with the PDF uncertainties has been presented in Ref. [18], which also takes into 
account the correlation between the scale and the PDFs (which are evaluated at a 
given factorization scale): one calculates the maximal/minimal cross sections with 
respect to the scale variation and apply on these cross sections the PDF-|-A''^Pa<j-|-A*'^acj 
uncertainty. This procedure has been in fact already proposed, together with other 
possibilities which give similar results, in Ref. [57] where top quark pair production at 
hadron colliders was discussed. An argument for the approach put forward in Ref. [57] 
is that unknown high order effects also enter in the PDFs determination and one 
needs to use the full information on the PDF uncertainties in the determination of the 
scale dependence. To this combined scale+PDF uncertainty, one can linearly add the 
scheme/EFT uncertainty to obtain the overall theoretical error. 

B) In fact, a much simpler way to combine the scale, EFT/scheme and the PDF uncer- 
tainties which, as mentioned before, should all be considered as pure theoretical errors, 
is to add them linearly. Doing so, one in general obtains uncertainties that are slightly 
larger than with the procedure A described above. 

C) A third possibility is to take the different predictions of the cross sections when eval- 
uated with the various sets of NNLO PDF sets as a measure of the PDF uncertainty. 
Since the maximal cross section for gg ^ H is obtained with the MSTW parametriza- 
tion and the minimal one with the ABKM set (see Fig. the total uncertainty on the 
production cross section will be simply the sum of the scale and scheme/EFT uncer- 
tainties when evaluated with the MSTW PDF set (for the -|- uncertainty) and when 
evaluated with the ABKM PDF set (for the - uncertainty) which maximizes/minimizes 
the cross section. 

^"in statistical language, the PDF uncertainties should be considered as having a flat prior, exactly like the 
scale uncertainty. A more elaborated discussion can be found in [14] where the recommended combination 
is that of the linear type. 
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Figure 7: The production cross section <j{gg — > H) at NNLO at tlie IHC with -/s = 7 TeV, 
including the total theoretical uncertainty band when all the individual uncertainties are combined 
using the three procedures A, B and C described in the text. In the inserts, the relative deviations 
from the central cross section value are shown. 

Using the three procedures A, B, C described above, one obtains the overall theoretical 
uncertainties on (T^^^'^{gg — t- H) at the 7 TeV IHC that are shown in Fig. |^ as a function 
of the Higgs mass. With the procedure A, the uncertainty amounts"*^^ to ~ —23%, +25% for 
Mh < 160 GeV, reduces to ^ -21%, +22% for Mh > 200 GeV to reach the value « ±22% 
at Mh ~ 500 GeV. With the procedure B, i.e. when the various sources of uncertainties 
are added linearly, one obtains approximately the same total uncertainty as above; they are 
a few percent (2 to 3%) smaller for the + uncertainty and higher for the - uncertainty. In 
the case of procedure C, as we use for the central prediction of a^^^°{gg -)> H) the MSTW 
PDF set, the + uncertainty is simply the sum of the scale and scheme uncertainties which 
add to less than 20% in the entire Higgs mass and are thus much smaller (in particular in 
the low Higgs mass range) than in procedures A and B, while for the - uncertainty, one 
has to add the difference between the MSTW and ABKM predictions. In this case, the 
overall uncertainty is comparable to what is obtained using the procedures A and B in 
the low Higgs mass range where the difference between MSTW and ABKM is of the order 
of 10% and, hence, the Hessian MSTW POF+a^ uncertainty. It becomes much larger 
at higher Mh values where the difference between the MSTW and ABKM predictions 
becomes significant: at Mh = 500 GeV, one obtains an uncertainty of ~ —35%, compared 
to ~ —23% using procedures A and B. 



^^Note that the overall uncertainty on a^^^'~'{gg — >■ H) at the IHC is significantly smaller than what has 
been obtained at the Tevatron using the procedure A for the combination , ~ —35%, +40%, as a result of 
the reduction of both the scale and the PDF uncertainties at — 7 TeV compared to y's = 1.96 TeV. 
The main reasons are that, at the IHC, the QCD /('-factors are smaller and the gluon density has a less 
uncertain behavior for not too heavy Higgs bosons. 
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Mh 


^ A B C 


115 


-lo or+9.1%+9.1%+7.3% +27.8% +25.6% +16.4% 
10.0d_^p 2%-8.8%-7.3% -24.6% -26.3% -26.8% 


120 


-If? D^+8.9%+9.1%+7.3% +27.5% +25.3% +16.2% 
lU.O'±_^p 2%-8.8%-7.3% -24.5% -26.2% -26.8% 


125 


1 r CI +8.8%+9.1%+7.2% +27.1% +25.1% +16.0% 
-^^•"^-^-9.9%-8.8%-7.2% -24.1% -25.8% -26.7% 


130 


14 q9+8.5%+9.1%+7.1% +26.7% +24.7% +15.6% 
J-'*-'5^_9.8%_8.8%-7.1% -24.0% -25.6% -26.7% 


135 


in 9f^+8.4%+9.1%+7.0% +26.4% +24.5% +15.4% 
■^"*-"'"-9.6%-8.8%-7.0% -23.7% -25.4% -26.6% 


140 


19 oi+8.3%+9.1%+7.0% +26.2% +24.3% +15.3% 
-^^■"^-^-9.5%-8.8%-7.0% -23.6% -25.2% -26.7% 


145 


1 1 4r+8.2%+9.1%+6.9% +26.0% +24.2% +15.1% 
-^^•^"'-9.5%-8.8%-6.9% -23.6% -25.2% -26.8% 


150 


-in R7+8-l%+9.1%+6.8% +25.8% +24% +14.9% 
-^"•"'-9.5%-8.8%-6.8% -23.5% -25.1% -26.9% 


155 


Q Q/l+7-9%+9.1%+6.6% +25.4% +23.6% +14.5% 
»-»'i_9.4%_8.8%-6.6% -23.2% -24.8% -26.8% 


160 


Q 91+7.8%+9.l%+5.9% +24.4% +22.8% +13.7% 
»-^-L_9.4%_8.8%-5.9% -22.7% -24.2% -26.4% 


165 


.7+7.7%+9.1%+4.9% +23.4% +21.7% +12.6% 
<'-*'-9.4%-8.8%-4.9% -21.7% -23.1% -25.5% 


170 


7 o7+7.7%+9.1%+4.2% +22.7% +21.1% +11.9% 
'■°'-9.4%-8.8%-4.2% -21.0% -22.5% -25.0% 


175 


7 or+7.6%+9.1%+3.7% +22.0% +20.4% +11.3% 
' •'^'-'-9.4%-8.9%-3.7% -20.5% -21.9% -24.6% 


180 


fi OG+7.5%+9.2%+3.1% +21.4% +19.8% +10.6% 
"•°"-9.3%-8.9%-3.1% -19.9% -21.3% -24.2% 


185 


(> /i9+7.4%+9.2%+3.0% +21.1% +19.6% +10.4% 
"•^^-9.3%-8.9%-3.0% -19.8% -21.2% -24.2% 


190 


ni +7.4%+9.2%+3.4% +21.5% +20.0% +10.8% 
"•"-^-9.3%-8.9%-3.4% -20.3% -21.7% -24.9% 


195 


r- f.r+7.4%+9.2%+3.6% +21.8% +20.3% +11.0% 
'^•"'-'-9.3%-8.9%-3.6% -20.5% -21.9% -25.2% 


200 


c q/i+7.3%+9.3%+3.7% +21.8% +20.3% +11.0% 
^•'^^-9.3%-9.0%-3.7% -20.6% -21.9% -25.4% 


210 


A n-| +7.2%+9.3%+3.7% +21.7% +20.2% +10.9% 
^•°-^-9.3%-9.0%-3.7% -20.6% -22.0% -25.8% 


220 


A QR+7.2%+9.3%+3.6% +21.6% +20.2% +10.9% 
^•"'"-9.2%-9.1%-3.6% -20.6% -21.9% -26.0% 


230 


n Q7+7.0%+9.4%+3.5% +21.4% +19.9% +10.5% 
O.y/ _g 2%-9.2%,-3.5% -20.5%: -21.8%, -26.2% 


240 


/.- + 7.()'X. + y..^.'X. + :S,:!'X +21.1% -1<).,S'/; +10.:!% 
'■'■'■'"'-9.2'X-I).2%,-:S.3% -20.5% -21.8%, -2().5%: 


250 


q .,7+6.9%+9.5%+3.1% +20.9% +19.5% +10.0% 
'^•'^'-9.2%-9.3%-3.1% -20.4% -21.6% -26.6% 


260 


q -1 -| +6.8%+9.6%+3.0% +20.6% +19.3% +9.7% 
0.1i_g 2%_g,4%_3,o% -20.3% -21.6% -26.8% 


270 


9 nQ+6.7%+9.7%+2.8% +20.5% +19.1% +9.4% 
Z.oa_;, 2%-9.5%,-2.8% -20.2% -21.4% -27.0% 


280 


9 7-,+6.8%,+9.8%,+2.6% +20.5% +19.2% +9.4% 
Z.(l_g 2%_9 5%_2.6% -20.1% -21.4% -27.2% 


290 


9 cr+6.8%+9.8%+2.4% +20.3% +19.1% +9.2% 
^•'JJ_9.1%_9.6%-2.4% -20.0% -21.1% -27.2% 


300 


9 49+6.7%+9.9%+2.3% +20.2% +18.9% +9.0% 
"'•^'^-9.1%-9.7%-2.3% -19.9% -21.1% -27.4% 


320 


9 9q+6.7%+10.1%+2.3% +20.2% +19.0% +9.0% 
^•^"'-9.2%-9.9%-2.3% -20.1% -21.4% -28.2% 


340 


9 -|q+6.9%+10.3%+3.0% +21.4% +20.2% +10.0% 
^•-^^-9.2%-10.1%-3.0% -21.1% -22.3% -29.6% 


360 


9 qi +7.0%+10.5%+4.1% +22.5% +21.6% +11.0% 
z.oi_g 2%-l0.3%-4.1% -22.3% -23.6% -31.3% 


380 


9 -io+6.3% + 10.7%+2.5% +20.4% +19.6% +8.8% 
^•-^°-9.1%-10.5%-2.5% -20.9% -22.2% -30.3% 


400 


-1 Qq+5.9%+11.0%+3.1% +20.7% +20.0% +9.0% 
-^•^"*-8.8%-10.7%-3.1% -21.5% -22.7% -31.3% 


450 


-1 97+5.0%+11.6%+4.0% +21.4% +20.5% +9.0% 
■L-^'_8.4%-11.3%-4.0% -22.6% -23.8%, -33.4% 


500 


n 7Q+4.4%+12.2%+4.5% +22.1% +21.1% +8.9%, 
'^•'^-8.1%-11.9%-4.5% -23.3% -24.5% -35.2% 


600 


q-|+3.7%+13.3%+6.6% +24.5% +23.7% +10.4% 
"•"'■^-7.7%-13.0%-6.6% -26.1% -27.3% -40.0% 



Table 2: The NNLO total Higgs production cross sections in the gg ^ H process at the IHC with 
■v/s = 7 TeV (in pb) for given Higgs mass values (in GeV) at a central scale jjLp = jJiR = \Mh- 
Shown also are the corresponding shifts due to the theoretical uncertainties from the various sources 
discussed (first from scale, then from PDF+A°'^p+"^Q!s at 90%CL and from EFT), as well as the 
total uncertainty when all errors are added using the procedures A, B and C described in the text. 
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In Table we summarize the results obtained in this section^^. The gg ^ H pro- 
duction cross section for values of the Higgs mass relevant at the IHC with y/s = 7 TeV 
are given together with the uncertainties from scale variations, the POF+A^^Pa^+^^'^Os 
uncertainty in the MSTW scheme and the uncertainty due the use of the EFT approach 
beyond the NLO approximation. The combined uncertainties obtained using the three 
procedures A, B and C proposed above are also given. 



2.4 The cross sections at ^ = 8 10 TeV and 14 TeV 

As already stated, intermediate center of mass energies between y/s = 7 and 10 TeV are 
being currently considered for the IHC in a very near future, with the possibility of collecting 
integrated luminosities that are significantly larger than 1 fb~^ [11]. We therefore also 
display the theoretical predictions for the Higgs production rates in the gg ^ H channel 
at these energies, following exactly the same steps as the ones discussed in section 2.1. 
The results for a^'^^'^{gg — H), again for a central scale = fip = \Mh and using the 
MSTW2008 set of NNLO PDFs, are displayed for ^ = 8,9 and 10 TeV in Fig. | as a 
function of Mh and in Table H for the relevant Higgs mass values . 



s = 8 TeV 
^ = 9 TeV 
v/s = 10 TeV 




150 200 300 

Mh [GeV] 



400 



500 



Figure 8: The production cross section (j{gg H) at NNLO as a function of Mh at the IHC with 
center of mass energies of ^ = 8,9 and 10 TeV. In the inserts, the relative increase compared to 
the cross section at -^s ~ 7 TeV are shown. 

We see that in the low Higgs mass range, Mh ~ 120 GeV, the cross section a{gg — )• H) 
is ^ 20%, 40% and 100% higher at, respectively, ^/s = 8, 9 and 10 TeV, compared to 
^/s = 7 TeV. At higher Higgs masses, Mh ^ 300 GeV, the increase of the cross section is 
slightly larg result of the reduction of the phase space at lower energies. 

As far as the theoretical uncertainties are concerned, we have verified that they differ 
only by a small amount, less than one to two percent, when increasing the center of mass 
energy from ^/s = 7 to y/s = 10 TeV. To a good approximation, one can therefore use the 



^An extended table for Higgs masses up to 1 TeV can be found in the (;(;-fusion section of Ref. [14]. 
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Mh 


8 TcV 


9 ToV 

'^gg-^H 


„W TcV 

'^gg-^H 


Mh 


8 TcV 

^gg^H 


9 ToV 

'^gg^H 


„W ToV 

'^gg^H 


115 


23.31 


28.63 


34.26 


210 


6.39 


8.15 


10.07 


120 


21.46 


26.42 


31.68 


220 


5.82 


7.44 


9.22 


125 


19.81 


24.44 


29.37 


230 


5.33 


6.84 


8.50 


130 


18.34 


22.68 


27.30 


240 


4.91 


6.32 


7.88 


135 


17.03 


21.10 


25.44 


250 


4.55 


5.88 


7.34 


140 


15.84 


19.67 


23.76 


260 


4.22 


5.47 


6.85 


145 


14.77 


18.38 


22.24 


270 


3.94 


5.13 


6.44 


150 


13.80 


17.21 


20.86 


280 


3.70 


4.83 


6.08 


155 


12.89 


16.10 


19.55 


290 


3.50 


4.58 


5.78 


160 


11.97 


14.98 


18.22 


300 


3.33 


4.37 


5.53 


165 


11.04 


13.85 


16.87 


320 


3.10 


4.09 


5.20 


170 


10.28 


12.92 


15.77 


340 


3.06 


4.07 


5.19 


175 


9.62 


12.11 


14.80 


360 


3.26 


4.35 


5.58 


180 


9.00 


11.36 


13.90 


380 


3.09 


4.15 


5.36 


185 


8.44 


10.66 


13.08 


400 


2.76 


3.73 


4.83 


190 


7.92 


10.03 


12.32 


450 


1.85 


2.53 


3.32 


195 


7.47 


9.47 


11.65 


500 


1.17 


1.62 


2.15 


200 


7.07 


8.99 


11.07 


600 


0.47 


0.68 


0.92 



Table 3: The cross sections in the a^^^^{gg -> H) at the LHC with = 8, 9, 10 TeV (in pb) for 
given Higgs mass values (in GeV) at a central scale /i_F = ^r, = using the MSTW PDF set. 

values of tlie scale, sclieme/EFT and PDF+A'^^Pas+A*'^as uncertainties, as well as the 
combined uncertainties, that are given in Table ^ for the ^/s = 7 TeV case. 

Finally, let us briefly summarize the expectations for the LHC with the design center 
of mass energy of -v/s = 14 TeV, which is expected to collect at least ~ 30 fb~^ of data that 
should allow to discover the Higgs boson in its entire mass range 115 GeV<Mj^ < 1 TeV. 



10 



The results for the gg ^ H cross section at ^/s = 14 TeV are summarized in Figs. 
11 and in Table ^, following exactly the same lines and procedure as at y/s = 7 TeV. 
We will simply make a few comments to highlight the main differences: 

- The scale uncertainty, estimated by varying /i^ and fip again in the domain ^Mjj < 
fJ-R = f^F < Mh, does not significantly change when comparing to the IHC case: one 
obtains a variation of ~ ±10% at low and ~ ±5% at high Higgs masses; see Fig. |^ (left). 

- The EFT/scheme uncertainty is almost exactly the same than at ^/s = 7 TeV, as its 
most important component enters as a multiplicative factor in the gg^H amplitude but 
is larger starting from Mh ^ 500 GeV; Fig. ^ (right). 

- For the PDF uncertainties, one can see in Fig. [l^ (left) that one still has significant 
differences when folding the partonic gg ^ H cross section with the gluon luminosities 
given by the four NNLO PDFs sets. Fig. |l| (right) displays the PDF, PDF+A^^Pa^ 
and the combined PDE+A'^^Pas+A^'^a^ uncertainties at the LHC. One obtains a slightly 
smaller PDF+a^ uncertainty than at the IHC, 1 to 2% , since lower x values are probed. 
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Mh [GeV] Mh [GeV] 

Figure 9: The uncertainty bands of cr^J^^ from scale variation witli n — 2 (left) and the total 
EFT uncertainty (right) at -^s = 14 TeV as a function of Mh- 




Mh [GeV] Mh [GeV] 

Figure 10: The PDF uncertainties in a^^^°{gg H) at the LHC with ^/s = 14 TeV as a function 
of Mh- Left: the central values when using the four NNLO PDFs and right: the 90% CL PDF, 
PDF+A'^'^Pas and PDF+A'^'^Pa^ + A'^a^ uncertainties in the IVISTW scheme. 



a(gg ^ H) [pb] 




Mh [GeV] 

Figure 11: The cross section a^^^'^{gg H) at the LHC with the uncertainty band when all 
theoretical uncertainties are added using our procedure A described in the text. 
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- For the overall uncertainty on cr^^^^{gg — >• H) shown in Fig. 11, where we have 
restricted ourselves to procedure A only, it is more or less the same at 14 TeV than at 7 
TeV in the low Higgs mass range, but is slightly smaller for heavier Higgs bosons. As an 
example, we obtain for Mh = 500 GeV a total w ±21% uncertainty at 14 TeV compared 
to ^ ±23% at 7 TeV. 

Table ^ which displays the cross sections together with the individual and overall 
theoretical uncertainties for the Higgs masses relevant at the LHC, summarizes our results. 



Mh 




Mh 


""gg^H A B 


115 


„„+g.4%-|-8.7%+7.8% +28.1% +25.8% 
oy-^'_12.2%-8.5%-7.8% -26.6% -28.5% 


240 


1 ^ qn+7.0%+8.0%+3.7% +20.0% +18.7% 
iO.^U_g 20/^_7 g%_3 7% -18.6% -19.6% 


120 


rjp._|_g.2%+8.6%+7.7% +27.8% +25.5% 
^^■'^^-ll,Q%-8A%-7.7% -26.1% -28.0% 


250 


-, ^ oo+6.9%+8.0%+3.5% +19.7% +18.4% 
i4.d5_g 2%_7 g%_3 5% -18.4% -19.5% 


125 


.^+9.o%+8.5%+7.6% +27.4% +25.1% 
^''■^^^11.6%-8A%-7.6% -25.7% -27.5% 


260 


„ (.„+6.8%+8.0%+3.3% +19.4% +18.1% 
-L'J-0^_8.l%-7.8%-3.3% -18.2% -19.2% 


130 


zlR /~iQ+8.9%+8.5%+7.5% +27.2% +25.0% 
^°-^^-llA%-8.3%-7.5% -25.5% -27.3% 


270 


w „ _„+6.7%+8.0%+3.1% +19.1% +17.9% 
-'-^■'^-8.1%-7.8%-3.1% -18.1% -19.1% 


135 


r,r.+8.7%+8.4%+7.5% +26.8% +24.6% 
40.U0_^^ -^,j,_g 2%_7.5% _25.1% -26.8% 


280 


, „ -. _+6.7%+8.0%+2.9% +18.9% +17.7% 
-'-^■-l-'_8.i%_7.8%-2.9% -17.8% -18.8% 


140 


49 or\+S-5%+8Ayo+7Ayo +26.4% +24.3% 
^^-JU^jq g%_g 2%_7.4% -24.7% -26.4% 


290 


^ „|.+6.6%+8.0%+2.8% +18.6% +17.4% 
ii.DO_g go/^_2 g% -17.6% -18.6% 


145 


on Qn+8.4%+8.3%+7.3% +26.0% +24.0% 
•jy-»^J_10.6%-8.1%-7.3% -24.4% -26.1% 


300 


^ fjrj+6.5%+8.0%+3.4% +19.2% +18.0% 
-L-l--^^_8.o%-7.8%-3.4% -18.4% -19.3% 


150 


„„ p.„+8.3%-|-8.3%+7.2% +25.8% +23.8% 
■°^-10.4%-8.1%-7.2% -24.1% -25.7% 


320 


„ „„+6.5%+8.1%+3.1% +18.8% +17.8% 
-'-'^■'^-8.0%-7.9%-3.1% -18.1% -19.1% 


155 


„j. oo+8.1%+8.3%+7.0% +25.3% +23.4% 
•JO-'J^_10.2%-8.1%-7.0% -23.7% -25.3% 


340 


, p. oq+6.5%+8.1%+2.8% +18.4% +17.4% 
-LU.OO_g Q0/^_7 g%_2.8% -17.8% -18.7% 


160 


„„ „g+8.0%+8.2%+6.3% +24.4% +22.6% 
•J'J-Uo_;^g po/^_g go/^_g _22.9% -24.4% 


360 


, _„+6.4%+8.1%+3.5% +19.0% +18.1% 
-L-L-' '_8.o%-8.0%-3.5% -18.6% -19.5% 


165 


„„ _,_,+7.9%+8.2%+5.3% +23.1% +21.4% 
'-9.8%-8.0%-5.3% -21.8% -23.1% 


380 


1 1 /lc+6.0%+8.2%+4.4% +19.5% +18.6% 
-'--'-■^"-7.7%-8.1%-4.4% -19.2% -20.2% 


170 


no oQ+7.8%+8.2%+4.5% +22.3% +20.6% 
^O-O»_g.7%_8.0%-4.5% -21.0% -22.3% 


400 


in 4/^+5. 6%+8.2%+5.0% +19.7% +18.9% 
-^'^■^"-7.4%-8.1%-5.0% -19.7% -20.6% 


175 


97 O/|+7.8%+8.2%+4.0% +21.6% +20.0% 
^' ■^^-9.5%-7.9%-4.0% -20.2% -21.4% 


450 


7 /,r,+5.0%+8.4%+6.0% +20.0% +19.3% 
'•*^-7.0%-8.3%-6.0% -20.4% -21.3% 


180 


or 7i+7.7%+8.2%+3.5% +20.9% +19.3% 
^^■'-^-9.4%-7.9%-3.5% -19.6% -20.8% 


500 


A Q7+4.6%+8.6%+6.4% +20.4% +19.7% 
^■^'-6.7%-8.6%-6.4% -20.9% -21.7% 


185 


94 9S+7-6%+8.1%+3.3% +20.6% +19.0% 
^^■^°-9.1%-7.9%-3.3% -19.2% -20.4% 


550 


q q9+4.3%+8.9%+7.4% +21.3% +20.7% 
"'■'^^-6.5%-8.8%-7.4% -21.9% -22.7% 


190 


99 Q7+7-6%+8.1%+3.8% +21.0% +19.5% 
^^■■^'-9.1%-7.9%-3.8% -19.6% -20.7% 


600 


9 94+4.1%+9.2%+8.3% +22.1% +21.6% 
^•^^-6.3%-9.1%-8.3% -22.9% -23.7% 


195 


91 S'?+''-5%+8-l%+4-0% +21.1% +19.6% 
^J--0'J_g.oo/„_7.9%_4.0% -19.7% -20.8% 


650 


1 cq+3.9%+9.5%+9.0% +22.7% +22.4% 
1.00_g 2%_g 4y^_g 0% -23.9% -24.6% 


200 


on oo+7.4%+8.1%+4.1% +21.0% +19.5% 
^'J-0'i-8.8%-7.9%-4:A% -19.6% -20.7% 


700 


1 nt:+3-8%+9.8%+9.6% +23.4% +23.2% 
-^■"^-6.1%-9.6%-9.6% -24.6% -25.3% 


210 


10 in+'i'-3%+8.1%+4.1% +20.9% +19.5% 
-^^■-^^-8.6%-7.8%-4.1% -19.4% -20.5% 


800 


n r::o+3.5%+10.4% + 10.5% +24.6% +24.4% 
"•'^^-6.0%-10.2%-10.5% -26.1% -26.7% 


220 


1 7 f?/l+7.2%+8.1%+4.0% +20.6% +19.3% 
-L' •O'±_8.4%_7.8%-4.0% -19.2% -20.3% 


900 


n 97+3.4%+11.0% + 11.3% +25.6% +25.7% 
"•^'-5.9%-10.7%-11.3% -27.5% -28.0% 


230 


1^; qn+7.1%+8.0%+3.8% +20.4% +19.0% 
l"-'JO_8.3%-7.8%-3.8% -18.9% -20.0% 


1000 


n 1 t;+3.3% + 11.5%+12.0% +26.4% +26.9% 
"•-^"-'-5.8%-11.3%-12.0% -28.8% -29.1% 



Table 4: The NNLO total Higgs production cross sections in the gg ^ H process at the LHC with 
^/s = 14 TeV (in pb) for given Higgs mass values (in GeV) at a central scale — fJ-R — \Mh- 
Shown also are the corresponding shifts due to the theoretical uncertainties from the various sources 
discussed (first from scale, then from PDF+A°''P+"^as at 90%CL and from EFT), as well as the 
total uncertainty when all errors are added using the procedures A and B described in the text. 
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3. The Higgs decay branching fractions 
3.1 The parametric uncertainties 

There is an additional source of theoretical errors which has not been considered in Ref. [18] 
and not included in the Tevatron and LHC experimental analyses [9, 12, 13]: the one due 
to the branching ratios in the Higgs boson decays^^. Indeed, contrary to the partial widths 
of Higgs decays into leptons and gauge bosons where only elcctroweak corrections appear 
at the lowest orders and are thus known precisely, the partial decays widths into quark 
pairs and gluons are plagued with various uncertainties due to strong interactions [23]. 

As it is well known, the QCD corrections to the H QQ partial widths, with Q = b 
and c, are extremely large; for reviews, see Refs. [2,6]. The bulk of these corrections can 
be absorbed into running quarks masses by defining an effective HQQ Yukawa interaction 
in which the pole quark masses Mq{Mq) are replaced by running quark masses in the MS 
scheme, mg, with ag evaluated at the scale of the pole mass fi = Mq [58]. One then evolves 
Wq from Mq upward to a renormalization scale ji, which should be chosen to be close to 
the Higgs mass, /x « Mh- Using as starting points the pole masses Mq, the values of the 
running b, c masses at the scale fi ~ 100 GeV are, respectively, ^1.5 and Ri 2 times smaller 
than the pole masses. The partial widths, r{H QQ) oc friQ (but where additional QCD 
corrections that are known up to three-loops need to be implemented [6,59]), are thus 
reduced by factors of ~ 2 and ~ 4 for Higgs decays into bb and cc, respectively. 

Nevertheless, residual uncertainties in the H — )• QQ partial decay widths will remain. 
These uncertainties are mainly due to three sources: i) the imperfect knowledge of the input 
bottom and charm quark masses, ii) the error on which migrates into uncertainties on 
the running quarks masses at the scale // and on the residual QCD corrections to the partial 
widths and, finally, Hi) the variation of the scale fi which enters as an argument in as, in 
the running quarks masses mg and in the residual QCD corrections. This last uncertainty, 
which is expected to account for the missing (not yet calculated) higher order contributions 
to the partial widths, turns out to be very small (as both the quark mass evolution and 
the residual corrections are known up to three loops) and can be safely neglected [60] . 

In the case of the decay H — t- gg, the leading order partial width, that is proportional 
to Q^, has to be corrected by large QCD corrections which nearly double the width [20, 61]. 
It is affected by uncertainties due to the error on and to the residual scale dependence 
which is about ~ 10% at NNLO [61]. However, as the branching ratio of this decay is 
of order a few percent at most in the entire Higgs mass range, these uncertainties do not 
affect the total Higgs width and thus, the dominant branching ratios, in a substantial way. 

We will assume that the uncertainties in the decays widths H — t- t'^t~ and H — t- VV* 
are negligible^^ as QCD effects will enter the partial widths only at higher orders, > O(a^). 

^^A discussion of QCD uncertainties in Higgs decays has also been made recently in Ref. [14]. 

^^This might not be entirely true in the case of the H — 5> WW and H — > ZZ decay channels near the 
2Mw and 2Mz thresholds where the regularization of the spurious spikes that appear when including one- 
loop electroweak corrections [62], might introduce some ambiguities. Furthermore, for heavy Higgs bosons 
Mh > 500 GeV, additional corrections to the partial H — >■ WW, ZZ widths from the Higgs self-coupling, 
\hhhoc.M% become very large [63]; these uncertainties mostly cancel in the branching ratios, though (but 
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In the case oi H ^ 77, the decay is mediated mainly by a loop and to a lesser extent 
by a top-quark loop; in the low to moderate Mh range where the decay is important, the 
uncertainty due to the latter contribution is small, less than 1%, as the QCD corrections 
are tiny and the large mt limit a good approximation [64]. 

3.2 Uncertainties on the Higgs branching ratios 

An analysis of the uncertainties affecting the Higgs branching ratios has been performed 
long ago in Ref. [23]. Here, we update this analysis in the following way: 
i) For the strong coupling constant as, we adopt for consistency reasons the same value 
that we use for determining the central predictions for the production cross section, 
that is, the MSTW value with its associated la error, as(M|) = 0.1171 ± 0.0014 at 
NNLO. 

a) As starting points for the 6, c masses we use the MS masses mg evaluated at the mass 
itself with values^s mb(mi,) = 4.19;[]:J^ GeV and mc{mc) = l-27tom GeV [54]. When 
we will quote separately the impact of these uncertainties on the Higgs branching ratio, 
we will assume the central value as{M^) = 0.1171 at NNLO for the strong coupling 
constant. 

in) For the H ^ gg channel we will consider only the parametric uncertainties that is, 
the error on Ug and quark masses (the latter being negligible). At NNLO, the scale 
uncertainty is of the order 10% but since the branching ratio is small it migrates into 
an error of less than 1%. In addition, the N^LO contribution would reduce this scale 
variation to the level of a few percent at most and this uncertainty can be safely 
neglected. 

Using an adapted version of the program HDECAY [66] which calculates the partial 
and total Higgs decay widths including the relevant QCD and electroweak higher order 
corrections, we have evaluated^® the uncertainties discussed above and added in quadrature 
the errors due the input masses mfe,mc and the coupling as{M'^). We obtain the branching 



ratios (BR) shown in Fig. 12 as a function of the Higgs mass, including the total uncertainty 
bands. In Table ^, the branching ratios for the main decay modes H ^ bb and H — )• 
WW are given for three selected values of the Higgs mass, Mh = 120, 135 and 150 GeV, 
together with the individual and total uncertainties. The uncertainties (in percentage) 



they will substantially affect the cross sections for Higgs production in the vector boson fusion processes 
which are related to the H — > WW/ZZ partial decay widths). 

^^While the central values of these masses are close to those obtained in Ref. [65] directly from QCD sum 
rules in a consistent 0{as) expansion, the corresponding uncertainties are larger. We thus conservatively 
adopt the PDG values [54], which could also account for possible theory errors; in any case the total 
uncertainties obtained when using the errors of Ref. [65] would nearly be equal to the uncertainty due to 
as error alone. Note that with these starting inputs and the value of as{Mz) given above, the central 
values of the 6, c pole masses are mi, — 4.71 GeV and Ale ~ 1-54 GeV with variation ranges of, respectively, 
4.64-4.90 GeV and 1.42-1.63 GeV. For the pole bottom-quark mass, the central value is rather close to 
the one adopted in the MSTW scheme for parameterizing the parton densities [38], mj — 4.75 GeV, and 
we expect that this difference will have no practical impact. 

^®We thank M. Spira and collaborators for pointing to us a numerical problem in the evaluation of — >■ gg 
with HDECAY in an earlier version of the paper. 
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on the branching ratios for the H — t- r+r^ and H — )• ZZ, 77 channels are the same as 
those affecting the H — >■ WW mode. In Table displayed are the branching ratios for 
the various decays as well as their corresponding total uncertainties for a selection of Higgs 
masses that are relevant for the IHC (and the Tevatron as well). 

The largest total errors are by far the ones that affect BR(i? — t- cc) which are of the 
order of 20%, while the errors on BR(i? — )• gg) are at the level of a few to 10% at most. In 
the case of the if — )• cc channel, it is mainly due to the uncertainty on the input charm- 
quark mass nic and, to a lesser extent, to the uncertainty on as ; their combination leads to 
a very strong variation of the charm-quark mass at the high scale /x, mc(/x) oc [as(/Lt)]^^/^^. 
For the H ^ gg channel the uncertainty is mainly due to the uncertainty on a^, and in both 
cases, the error on the input 6-quark mass leads to a 2 to 5% uncertainty. Nevertheless, 
since the branching ratios for these two decays are small, at most a few percent in the 
Higgs mass range of interest, the associated uncertainties will affect the Higgs total width, 
and hence the branching ratios for the other decay channels, in a less significant way. 



channel 



H ^bb 



H WW 



M 



H 



120 
135 
150 



120 
135 
150 



BR(%) 



65.40 
41.01 
16.07 



14.06 
39.86 
69.45 



ATOc Ato;, AqIs 



+0.7% 
-0.6% 

+0.4% 
-0.4% 

+0.2% 
-0.1% 



^3.4% 
-1.2% 

^6.0% 
-2.1% 

H8.7% 
-3.0% 



+0.7% 
-0.8% 

+ 1.3% 
-1.3% 

+ 1.9% 
-1.9% 



+0.7% 
-0.6% 

+0.4% 
-0.4% 

+0.2% 
-0.1% 



h2.3% 
-6.3% 

H.4% 
-4.0% 

FO.5% 
-1.6% 



+ 1.4% 
-1.4% 

+0.9% 
-0.9% 

+0.3% 
-0.3% 



ABR 



+3.6% 
-1.6% 

+6.2% 
-2.5% 

+8.9% 
-3.6% 



+2.8% 
-6.5% 

+1.7% 
-4.1% 

+0.7% 
-1.6% 



Table 5: The Higgs decay branching ratio BR(i7 — > bb) and BR(iJ -> WW) (in %) for given Higgs 
mass values (in GeV) with the corresponding uncertainties from the various sources discussed in 
the text; the total uncertainties (adding the individual ones in quadrature) are also shown. 

The uncertainty on the H ^ bb decay channel (which is important at the Tevatron) 
depends strongly on the considered Higgs mass range. For Mjj < 120 GeV where the 
branching ratio is by far the dominant one, BR(if — )• bb) > 65%, the uncertainty is less 
than a few percent -1%, -^2% at Mh = 100 GeV and « -2%, +4% at Mh = 120 GeV). 
This is mainly due to the fact that as the channel is dominant, the bb partial width is the 
major component of the total width and the errors partly cancel in the branching ratio. 
There is, however, a residual error stemming mainly from the input bottom mass and to a 
lesser extent, Og and the charm mass. 

In the mass range M//~ 120-150 GeV, the partial widths for H^bb and WW decays 
have the same magnitude and the two branching ratios have larger uncertainties. The 
uncertainty on BR{H^bb) increases from a few percent at MjyRil20 GeV to w^%,-|-9% 
at Mj|/~150 GeV. At the same time, the uncertainty on BIl{H ^WW), which is the same 
as the one on BR{H-^ZZ*) and BR{H-^jj), drops from ^-7%,+3% at Mh~120 GeV 
to ~— 2%, 4-1% at M//~150 GeV. Here, nic has little impact as BR(i7— t-cc) is small. 

For higher Higgs masses, Mh > 160 GeV, the Higgs width is mainly controlled by the 
H — )• WW decay (and H — )• ZZ for Mjj > 2Mz) and the uncertainty on the branching 
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Mh [GeV] 

Figure 12: The Higgs decay branching ratios as a function of Mh including the total uncertainty 
bands from the la errors on the input quark masses and the coupling as (the individual errors have 
been added in quadrature). 

ratio BR(iJ — )• WW) drops to a level below 1%. The uncertainty on BR{H — )• bb) stays 
relatively large, ~ —4%, +11%, but the branching ratio itself is too small to be relevant. 

For Higgs masses beyond Mh > 350 GeV, the decay channel H ^ ti opens up and has 
a branching ratio of ^ 20% at Mh ~ 500 GeV dropping to less than 10% for Mh < 400 
GeV or Mh > 700 GeV. It will be affected by the uncertainties on mt, Os as well as by some 
electroweak contributions from the Higgs self-couplings. This leads to a total uncertainty 
that is estimated to be below the 5% level [14], which translates to an uncertainty of less 
than 1% in the branching ratios for the important decays H — )• WW, ZZ. 

Thus, the uncertainties on the important Higgs branching ratios BR(// — )• WW, ZZ, 77) 
as well as BR(// — t- bb) can be significant in the intermediate mass range Mh ~ 120-150 
GeV where the bb and WW decays are competing with each other. These uncertainties 
should therefore be included in the experimental analyses. 

3.3 Combination with the cross section uncertainty 

In the previous section we have discussed the overall uncertainty on the Higgs production 
cross section. To this overall uncertainty, one has to add the uncertainty on the Higgs 
branching ratios obtained when considering a specific Higgs decay channel. This combi- 
nation^^ is in general rather complicated to perform when several production and decay 
channels are involved as is the case for the Tevatron where both qq — >• HV and gg ^ H 
production and H ^ bb and H — )• WW decays channel have to be considered, the main 
reason being that the rates for the H ^ bb and H — )• WW decays are anti-correlated, the 
sum of all branching ratios being equal to unity. At the IHC, the situation is much simpler 

^^Here, we will not address two (connected) issues which become important for heavy Higgs bosons: 
i) the effect of the total Higgs decay width which becomes large and ii) the interference between the 
gg ^ H ^ VV signal and the gg — >■ VV background. A study of these two effects is in progress [67]. 
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Mh 


bb A'°* 


CC A*°' 




WW zz TT 77 A'°* 


100 


79.48 +^.2% 


q on +20.9% 
O.yu _22.7% 


fi 04 +5.6% 


1.09 0.11 8.22 0.16 


105 


77 69 +2.4% 
u.oy -1,1% 


q 01 +20.9% 
0.01 -22.7% 


7 50 +5-^^° 


2.39 0.21 8.11 0.18 


110 


74 96 +2-6% 

11. M -1.2% 


q 00 +20.8% 
0.00 _22.7% 


7 07 +5.3% 
(.y< -s,Y%^ 


4.75 0.43 7.89 0.19 ^7^% 


115 


70 95 +3 1% 


q 40 +20.8% 
'J-'io -22.7% 


97 +5.1% 
-8.3% 


8.54 0.86 7.53 0.21 +^°^" 


120 


65.40 


q 91 +20.7% 
"^■^-^ -22.8% 


n q/l +4.9% 


14.06 1.57 7.00 0.22 +"^" 


125 


58.32 +t:8l 


9 Of? +20.7% 
Z.OO _22.8% 


Q 1 1 +4.6% 
^■-^^ -7.1% 


21.35 2.62 6.29 0.23 +?S 

— . y /o 


130 


50.00 tl-'il 


9 +20.6% 
^..10 _22.9% 


7 ^f! +4.2% 
'•00 -6.4% 


30.13 3.95 5.43 0.22 tli% 


135 


41 01 


9 rv-i +20.5% 
^•"^ -23.0% 


79 +3.9% 
"■'^ -5.5% 


39.86 5.42 4.49 0.21 ^I'lyo 


140 


on no +7.2% 
OZ.UZ _2.8% 


-, +20.5% 

1.0 ( -23.1% 


r 00 +3.6% 
O.DS 


49.92 6.82 3.53 0.19 ^1^% 


145 


23.60 l^;!^: 


1 -1^ +20.6% 
l^-l" -23.1% 


4 fil +3-2% 


59.82 7.86 2.62 0.17 ^2°% 


150 


16.07 tin 


n 70 +20.6% 
U-'y -23.2% 


q qi +3.0% 


69.45 8.21 1.79 0.14 tii% 


155 


q 47 +9-7% 


n 46 +20.7% 
-23.3% 


nq +2-8% 
z.uy _3po/^ 


79.28 7.33 1.06 0.10 tii% 


160 


o ro +10.4% 

0.00 -3,g%„ 


n 1 +20.7% 
-23.4% 


n 04 +2.7% 
0.84 _2^^<^^ 


90.63 4.19 0.40 0.05 ^04% 


165 


1 90 +10.7% 
-L-^-S _4.o% 


06 +20-8% 

U.UO _23.4% 


SO +2.6% 


95.97 2.22 0.14 0.02 ^01% 


170 


n 09 +10.7% 


04 +20-8% 
^■'-'^ -23.4% 


n 91 +2.6% 
^■^-^ -2.5% 


96.43 2.35 0.09 0.02 +°;°| 


175 


Q (50 +10.7% 

U.DO -4,1%, 


r, no +20.8% 
^■'-'"^ -23.5% 


r, 1 7 +2.5% 

-2.6% 


95.84 3.21 0.07 0.01 to°i% 
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r> (^9 +20.8% 
^■'-'^ -23.5% 


n 1 f\ +2.5% 


93.28 5.94 0.06 0.01 +°°| 
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U-^U -4,0%^ 


n no +20.8% 
^■^^ -23.4% 


n 1 9 +2.6% 
U-i^ -2.4% 
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190 


n 09 +10.7% 

U.O^ -4.1% 


n no +20.8% 

^■'-'^ -23.4% 


10 +2.5% 


78.72 20.77 0.04 0.01 +°;°| 


195 


n 90 +10.8% 


n ni +20.8% 

^■^1 -23.4% 


0.10 tin 


75.89 23.66 0.03 0.01 
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r, 9r +10.8% 

U./O _4.i% 


n ni +20.8% 

^■'-'-^ -23.5% 


0.09 


74.26 25.34 0.03 0.01 +°;°1 



Table 6: The Higgs decay branching ratios (in % ) for given Higgs mass values (in GcV) with the 
corresponding total uncertainties from the various sources discussed in the text. 

as, in practice, only the gg ^ H production channel is to be considered and only the decays 
H — )• WW, ZZ and to a lesser extent — ?• 77 are relevant. Since we have considered only 
the dominant QCD uncertainties, these decays are affected by the same uncertainties as 
discussed previously (see the last column of Table ^) . 

Nevertheless, combining the uncertainty in the production cross section and the uncer- 
tainty in the Higgs decay branching ratios is still not obvious. Indeed, while the uncertain- 
ties on the Higgs production cross section are purely theoretical, one should in principle 
consider the uncertainties on the Higgs branching ratios as experimental errors, since they 
are mainly due to the "experimental" determination of the quark masses and (and we 
have not included the corresponding theoretical errors). In addition, there is one param- 
eter which is common in the calculation of the branching ratios and the production cross 
sections: the coupling a^. The uncertainty on as will affect at the same time cr(gg — t- H) 
and BR(i7 — )• VV) and it occurs that, in both the Higgs production and in Higgs decays. 
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the minimal (maximal) values are obtained with the minimal (maximal) value of when 
the error A'^^^Os is included (we have assumed that no theoretical uncertainty on Og is 
present in the decays). 

Here, we will thus adopt the simple attitude of adding linearly the errors on the Higgs 
branching ratios to the theoretical uncertainty on the production cross section. In the case 
of the procedure A for the combination of the cross sections uncertainties, for consistency 
reasons, one should also use the 90% CL error on the branching ratios due to A^^^a^ 
together with the 90% CL POF+A^'^Pas uncertainties in the cross section. In this case, 
one obtains the combined uncertainty of a^^^'~'{gg H) x BR{H — )• VV*) at the IHC 
that is displayed in Fig. 13 in the case of the H — )• WW decay as a function of Mh- Besides 
the uncertainty on the cross section which is shown by the dashed lines, we display the 
effect of adding the error on the H — )• WW branching ratio as shown by the full lines which 
slightly increase the overall uncertainty. This is only visible in the region Mh ~ 100-150 
GeV where the errors on BR(// — )• WW), dominated in practice by the errors in the decay 
channel H — t- bb, is significant. Above the Mh ^ 2Mw threshold, the branching ratio 
uncertainty is small, below ~ 1%, and can be safely neglected^^. Note, however, that 
for large Mh values, the total Higgs width becomes significant and should be taken into 
account in both the production and in the decay [14]. 

In the case of the H — t- ZZ, 77 (and rr) decays, since the errors on the branching 
ratios are the same as those affecting BR(if — )• WW), the overall uncertainties in a x BR 



can also be seen from Fig. |13|: only the normalization of the branching ratio is different 
and can be obtained for a specific decay mode from Table |^. 



c7(gg ^ H) X BR(H ^ WW*) [pb] 

^ = 7 TeV 




120 130 



140 150 160 170 
Mh [GeV] 



190 200 



Figure 13: The production cross section times branching ratio for the process — > — ?> WW at 
the IHC with -^s = 7 TeV, including the total theoretical uncertainty band when the errors in the 
decay branching ratios are taken into account. In the inserts the relative deviations are shown. 



For instance, the branching ratio uncertainties have no eftect in the Higgs mass range Mh ~ 150-180 
GeV to which both the Tevatron and the IHC are most sensitive. However, at the Tevatron these errors 
should be taken into account for AIh < 150 GeV in both the H ^ bb and WW channels. 
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4. The MSSM neutral Higgs particles 



4.1 The MSSM Higgs sector 

In supersymmetric extensions of the SM, at least two Higgs doublet fields are required for 
a consistent electroweak symmetry breaking and in the minimal model, the MSSM, the 
Higgs sector is extended to contain five Higgs bosons: two CP-even h and H, a CP-odd A 
and two charged Higgs particles [4-7]. Besides the four masses, the properties of the 
MSSM Higgs sector are in principle determined by two more parameters: a mixing angle 
a in the neutral CP-even sector and the ratio of the vacuum expectation values of the two 
Higgs fields tan^; the value of the latter is expected to lie in the range 1 < tan^ < mt/nn,. 
In fact, in the MSSM, the scalar potential does not involve arbitrary self couplings as is 
the case with the SM, but involves only the gauge couplings and, as a result, only two free 
parameters are needed at tree-level to describe the Higgs sector; these are usually taken to 
be the mass Ma of the pseudoscalar A boson and tan p. In addition, while the masses of 
the heavy neutral and charged H, A, particles are expected to range from Mz to the 
SUSY breaking scale Ms = 0{\ TeV), the mass of the lightest Higgs boson h is bounded 
from above, M/j < Mz at the tree level. 

It is well known that this relation is altered by the large radiative corrections which 
enter in the Higgs sector, the leading part of which grow as the fourth power of mt and 
logarithmically with the SUSY scale or common squark mass M5; the mixing (or trilinear 
coupling) in the stop sector At and, to a lesser extent a Higgs mass parameter /x, play an 
important role. When higher order corrections are included, the upper bound on the mass 
of the lightest h boson is shifted from the tree level value Mz to the value MJ^^^ ~ 110-135 
GeV [7] in the maximal mixing scenario where Xf = At — /u/tan /3 \/^Ms. The maximal 
mixing scenario, together with the no-mixing scenario with Xt ^ which leads to a lower 
Mf^^^ value, are often used as benchmarks to describe MSSM Higgs phenomenology [68]. 

For a heavy enough A boson. Ma S> M^, the h boson reaches its maximal mass value 
M/i c:i Mjf^^ and has SM-like couplings to fermions and gauge bosons. In this decoupling 
regime which, for large tan/3 values occurs already for Ma ^ Mf^^^^, the three other Higgs 
particles are almost degenerate in mass, Mh ^ Ma^ M^. The couplings of the CP- 
even H boson (as well as those of the charged bosons) become similar to that of the 
pseudoscalar A boson: no tree-level couplings to the gauge bosons and couplings to isospin 
down (up) type fermions that are (inversely) proportional to tan /3. In particular, for high 
tan /3 values, tan /3 > 10, the H, A Yukawa couplings to bottom quarks and tau leptons are 
strongly enhanced, while those to top quarks are strongly suppressed. 

For a light pseudoscalar boson, Ma ^ M™^^ at high tan f3, the roles of the CP-even h 
and H states are reversed: it is the H boson which has a mass Mh ~ MJ^^ and SM-like 
couplings, while the h particle behaves exactly like the pseudoscalar A state, i.e. M^ 2± Ma, 
no couplings to gauge bosons and enhanced (suppressed) couplings to b, r (t) states. We are 
thus always in a scenario where one has a SM-like Higgs boson H^m = h (H) and two CP- 
odd like Higgs particles ^ = A ot H (h) when we are in the decoupling (anti-decoupling) 
regimc"*^^ which, for values tan/3 > 10, occurs already for Ma > M™^^ {Ma ^ MJ^^^^). 
'^^Note that there is an intermediate scenario for tan /3 > 10: the intense coupling regime [69] in which 
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Let us now discuss the Higgs Yukawa couplings to bottom quarks, which play a major 
role in the analysis that will be presented in this section. As mentioned earlier, in almost the 
entire parameter space for large enough tan (3 values, the couplings of one of the CP-even 
Higgs particles are SM-like, thus leading to the phenomenology described in the previous 
sections for the SM Higgs boson, while the couplings of the other CP-even particle are 
the same as those of the pseudoscalar A boson, on which we will focus in the rest of our 
discussion. At the tree level, this coupling is given in terms of the 6-quark mass, the SM 
vacuum expectation value v and tan /?, by 

\/2m6 tan/3>l \/2m6 a rf, A U fU\ ^AT\ 

A$f,fe = ^ — y tan(j , ^ = A,H [h] (4.1) 

V cos p V 

First of all, in the MSSM, one usually uses the modified dimensional reduction DR scheme 
which, contrary to the MS scheme, preserves Supersymmetry. In the case of the 6-quark 
mass, the relation between the DR and MS running masses at a given scale n reads [70] 



mD^(/.) = mr(^) 



3 TT vr^ 



(4.2) 



where the strong coupling constant as is also evaluated at the scale fi and additional 
but small electroweak contributions are present. Since the difference between the quark 
masses in the two schemes is not very large, Am^/m;, ~ 1%, to be compared with an 
"experimental" error on r?T5(mt) of the order of a few percent, it is common practice to 
neglect this difference, at least in unconstrained SUSY models where one does not have to 
evolve the parameters up to very high scales. We will thus adopt this approximation at 
least when we quote the central values of the cross sections (see later however, when the 
theoretical uncertainties will be considered). 



A second important fact is that the Yukawa coupling eq. (4J), besides the standard 
QCD corrections, receives additional one-loop vertex corrections which involve squarks and 
gluinos in the loops and which can be very large as they grow with tan /3 [71] , 

Afc — ^/UTTig tan/3/max(m^, ,mj ) (4.3) 

with m^^ ^ and rrig the sbottom and gluino masses; large corrections, proportional to ^At 
also appear when the electroweak stop and chargino loops are included. This correction 
is very important for large values of tan /3 and /U and significantly increase or decrease 
(depending on the sign of ^) the Yukawa coupling. It can be properly taken into account 
using a resumation procedure where one writes the modified Yukawa coupling, normalized 
to the coupling of the SM Higgs boson, as [70] 



tan /3 



1 + A. 



1 



tan2/3 



tan/3>l tan;3 



1 + A. 



(4.4) 



the three neutral Higgs bosons have comparable masses, Mh ~ Mh ~ Ma ~ M'^'^^ . The couplings of the 
three states to isospin down type fermions are enhanced and the squares of the CP-even Higgs couplings 
approximately add to the square of the CP-odd Higgs coupling. This property, together with the fact that 
Mh ~ Mh, makes that for our main purpose in this paper, one has the same phenomenology as in the 
decoupling or anti-decoupling regimes. 
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It has been shown that in this case all terms of 0{as tan f3) in perturbation theory are 
resummed [70]. The inclusion of the NNLO corrections in Af, allows, in addition, to have 
a very precise description of the coupling, with a scale variation at the percent level [72] . 

The resummation procedure which allows to include properly the SUSY radiative cor- 
rections in the ^bb Yukawa couplings is implemented in the major codes which calculate 
the MSSM Higgs spectra such as HDECAY [66], FeynHiggs [73] or the RGE codes that cal- 
culate the MSSM particle spectrum such as SuSpect [74]. Since in our analysis we will 
be mainly concerned with the standard QCD effects and we would like to have a model 
independent approach (see next subsection), we will assume that these SUSY QCD and 
electroweak corrections are handled by the previously mentioned programs that provide 
the radiatively corrected Higgs couplings. Nevertheless, we will see later that in the main 
channels discussed in this paper, the major part of these corrections cancels out. 

Before closing this survey of the MSSM Higgs sector, we note that there are experimen- 
tal constraints on the Higgs masses, which mainly come from the negative LEP2 searches [8]. 
In the decoupling limit where the h boson is SM-like, the constraint ^114 GeV holds 
and rules out tan /3 values smaller than ~ 3. Outside this regime, the most conservative and 
model independent limits on the Higgs masses from LEP2 searches are Mh, Ma ^ Mz [8]. 
Some exclusion limits in the [MA,tan/3] parameter space are also available from Tevatron 
searches [10] and rule out extremely large tan (3 and low Ma values. 

4.2 The production cross sections at the IHC 

In the MSSM, the dominant production processes for the CP-even neutral h and H bosons 
are essentially the same as those for the SM Higgs particle discussed in section 2.1. In fact, 
for the lighter h boson in the decoupling and the heavier H boson in the anti-decoupling 
regimes, the cross sections are almost exactly the same as for the SM Higgs particle with a 
mass ~ M™^^ shown in Fig. |l| with the numerical values of Table |l|. The only difference is 
that some SUSY particles might contribute to the gg fusion processes and, eventually, to 
the SUSY radiative corrections in the other processes. However, for heavy enough squarks 
and gluinos, as is the case for the benchmark scenarios that are in general discussed and 
which have Ms 3> M™*^^, these contributions are rather small as the SUSY states decouple 
from the amplitudes (their couplings to the Higgs bosons are not proportional to their 
masses). These SUSY contributions, including even the known QCD corrections [75], can 
be nevertheless evaluated exactly using the program HIGLU for instance. 

In the case of the CP-odd like particles, the situation is completely different. Because 
of CP invariance which forbids A couplings to gauge bosons at tree-level, the pseudoscalar 
A boson cannot be produced in the Higgs-strahlung and vector boson fusion processes; only 
the gg ^ A fusion as well as associated production with heavy quark pairs, qq, gg — )■ QQA, 
will be in practice relevant (additional processes, such as associated production of CP-even 
and CP-odd Higgs particles, have too small cross sections). This will therefore be also the 
case of the CP-even H and h particles in, respectively, the decoupling and anti-decoupling 
scenario. In addition, if one concentrates on the high tan /3 regime (that is the only relevant 
one in this context, both at the Tevatron and the IHC), the 6-quark will play the major 
role as its couplings to the CP-odd like bosons are enhanced. 
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Therefore, for ^ = A or H{h), one first has to take into account in the gg ^ ^ 
processes the &-quark loop which provides the dominant contribution in the MSSM, as the 
one of the top-quark loop is suppressed, X^tt oc 1/ tan^. Moreover, in associated Higgs 
production with heavy quarks, the rates for the pp tt^ processes are also suppressed 
and, instead, associated Higgs production with bb final states must be considered. In fact, 
the pp — >■ 66$ production processes become the dominant ones in the MSSM. 

In the gg ^ ^ processes with only the 6-quark loop included, as the b qiiark mass is 
very small compared to the Higgs masses, chiral symmetry approximately holds and the 
cross sections are approximately the same for the CP-even H (h) and CP-odd A bosons^°. 
The QCD corrections are known only to NLO for which the exact calculation with finite 
loop quark masses is available [20]. Contrary to the SM case, they increase only moderately 
the production cross sections. The calculation of the higher order corrections that have 
been made available for the SM Higgs boson, the NNLO QCD corrections (performed in 
the infinite quark mass limit) and the NLO electroweak corrections (the dominant part of 
which arises because of the large Higgs-tt Yukawa coupling) do not apply here and will 
be thus ignored. However, in the case of h and H production, to approach properly the 
decoupling and anti-decoupling limits and to reproduce the SM Higgs results discussed in 
section 2 when these higher order corrections are included, we will adopt the central value 
fiQ = T^Mq, for the rcnormalization and factorization scales, i.e. as in the SM case. 

In the case of the pp — )■ 66"I> processes, the NLO QCD corrections have been calculated 
in Ref. [76] and turn out to be rather large, in contrast to pp — > tt+Higgs production. 
Because of the small mf, value, the cross sections develop large logarithms log(Q^/m^) 
with the scale Q being typically of the order of the factorization scale, hf ^ ^ mi,. 
These logarithms can be resummed via the Altarelli-Parisi equations by considering the 6- 
quark as a massless parton and using heavy quark distribution functions at a scale jip ^ Q 
in a five active flavor scheme. In this scheme, the inclusive process where one does not 
require to observe the 6 quarks is simply the 2^1 process 66 ^ $ at leading order [24]. 
If one requires the observation of a high-pr Anal 6-quark, one has to consider its NLO 
corrections [25] and in particular the 2 — 2 process gb $6, which indeed generates the 
Pt of the 6-quark. Requiring the observation of two 6 quarks in the final state, one has 
to consider the 2^3 process gg — )■ 66$ discussed above, which is the leading mechanism 
at NNLO [26]. Thus, instead of qq^gg — )• 66$, we will consider the process bb ^ ^ for 
which the cross section is known up to NNLO in QCD [25, 26], with corrections that are of 
moderate size if i) the bottom quark mass in the Yukawa coupling is defined at the scale 
M$ to absorb large logarithms log(/x^/m^) and ii) if the factorization scale, that we will 
set here equal to the renormalization scale, is chosen to be small, = fxji = = |M$. 

In order to evaluate the Higgs production cross sections at IHC energies in these two 
main processes, gg ^ ^ and 66 — > $, we will proceed as follows. We only evaluate the cross 
sections for the pseudoscalar A boson: in the gg ^ A process at NLO using the program 
HIGLU [37] with a central scale fJ-R = fip = t^o = (and where only the bottom quark 

^°This is only true if the SUSY particle loop contributions are not included. In the case of the CP-even 

particles, their relative contribution are suppressed at high enough tan /3 and we will ignore them here. In 
the case of the A boson, the SUSY contributions appear only at two-loops and they can be safely neglected. 
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loop contribution is included by setting Xau = 0), and in the bb ^ A process up to NNLO 
using the program bbhSnnlo [26,77] with a central scale IJ^r = jJip = l^o = \Ma- In 
both cases, we work in the MS scheme for the renormalization of the bottom-quark mass. 
However, while rnhijni,) is used in the gg fusion process, rnh{fJ.R) is adopted in the bb fusion 
channel. In both processes, we assume the Abb coupling to be SM-like, that is, we will 
not include the tan /3 term and the SUSY corrections to Am;,. To obtain the true cross 
sections for A production, one will have therefore to rescale the numbers that we provide 



by the factor of eq. {iA) squared. Furthermore, to obtain the cross section for both CP- 
even and CP-odd Higgs production, an additional factor of two has to be included. As a 
consequence of chiral symmetry for M$ ^ and since the H (h) masses and couplings 
are very close to those of A, this turns out to be an excellent approximation. 

Within this set-up, the best values of the cross sections for gg ^ ^ and 66 — )• <1> are 
shown in Fig. ^ as a function of the Higgs mass M$ when the MSTW sets of (NLO for 
the former and NNLO for the latter process) PDFs are used to parametrize the gluon and 
bottom-quark densities. Center of mass energies in a range between = 7 TeV and 14 
TeV relevant for the LHC are considered. One can first notice that the cross sections for 
gg ^ ^ and 66 — t- <I> are comparable at a given energy and they significantly increase with 
increasing center of mass energies or decreasing Higgs mass. If, for example, the value 
tan/3 = 10 is adopted, the numbers in Fig. |l4| have to be multiplied by a factor ~ 200 to 
obtain the true cross section for both A and H{h) production. For low to moderate M$ 
values, the expected event rates are thus simply huge at the IHC, despite of the relatively 
low luminosities that are expected. This explains why the chances for observing a Higgs 
boson at the IHC are much higher in the MSSM than in the SM, as in the former case the 
production rates can be two to three orders of magnitude larger. 
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Figure 14: The production cross sections in the processes — ?> <& (left) at a central scale fi^ — 
IJ'F = fJ-o = and 66 — $ (right) at a central scale fiji — jip — — |Af$ as a function of M$ 

for several center of mass energies relevant for the LHC. The MSTW sets of PDFs at the required 
perturbative order have been used. Only the cross section in the pseudoscalar A case but with a 
SM-like Yukawa coupling is included. 
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4.3 The theoretical uncertainties 

Left to evaluate are then the theoretical uncertainties on the production cross sections and, 
for this purpose, we will follow very closely the procedure (A) discussed in section 2.2 for 
the SM Higgs boson, at least for the scale and PDF+Os uncertainties. For the numerical 
analysis, we will only consider the case of the IHC at ^/s = 7 TeV: the uncertainties at 
center of mass energies slightly above this value, y/s = 8-10 TeV, and even for the LHC 
energy = 14 TeV are expected to be comparable. The procedure and the main results 
that we obtain are briefly summarized in the following^^. 

In the case of the gg ^ ^ process, as for the SM Higgs boson, the scale uncertainty is 
evaluated by allowing for a variation of the renormalization and factorization scales within 
a factor of two around the central scale, ^/io < f^R^l^F < 2^o with jiQ = ^M$. For the 
66 — )• ^' case, we will extend the domain of scale variation to a factor of three around 
the central scale /Uq = l^f^, i/io < fJ'R, fJ'F < S/xq. The reason for this non democratic 
treatment is that, first, in the gg ^ ^ process, there is another uncertainty due to the 
scheme dependence in the renormalization of the bottom-quark which will add to the 
scale uncertainty, as will be discussed later. A second reason is that it is well known that 
when the same final states are considered, the cross sections in the 66 —t- $ process in 
the five-fiavor scheme and in the qq,gg — s- 66<I> channel in the four-fiavor scheme differ 
significantly [78] and only by allowing a wider domain for scale variation and, hence, a 
larger scale uncertainty that the two results become consistent with each other. In our 
final numbers for the uncertainties, we will assume a variation of ^fi and fip that is taken 
to be independent but with the additional restriction 1/k < fiR/fip ^ ^ imposed. However, 
to illustrate the much larger scale uncertainty that is possible in the 66 — t- <I> case, we will 
also show results when this constraint is relaxed, in much the same way as in Ref. [14]. 

The results for a{gg — )• $) and fT(66 — )• <I>) at the IHC, for scale variations in the 



domains fJ-o/n < fiR,HF < nfiQ with k = 2 and 3 are shown in Fig. 15 as a function of 
M$. One can see that the scale variation is moderate for gg ^ ^ with k = 2, despite of 
the fact that the process is known only at NLO, leading to an uncertainty of order ±10% 
in the entire Higgs mass range. Extending the variation domain to k = 3 will increase 
the uncertainty by another ~ 10%. As in the SM Higgs case, the maximal and minimal 
values of the cross sections are approximately obtained for /i/j ^ fip and thus, varying 
independently the two scale does not affect the uncertainty. For k = 2 when the restriction 
^ < IJ'r/ ^J-F < 2 is imposed, the uncertainty is also small in the case of 66 — )• $, « ±10% 
at low Higgs masses and less for higher masses; this was to be expected as the process 
is evaluated at NNLO. However, the extension to «; = 3, while keeping the restriction 
5 < /^fl/AiF < 3, will significantly increase the uncertainty: one would have —18%, ±24% 
at M<t, = 100 GeV and -13%, ±6% at M$ = 200 GeV. If the restriction on hr/iif is 
ignored the uncertainty blows up, especially for very low Higgs masses: one would have a 
variation of -45%, ±25% at M$ = 100 GeV, as also noticed in Ref. [14]. 

^^Note that since in our analysis we are not considering the SUSY particle contributions and focus only 
on the standard QCD effects, additional uncertainties from the SUSY sector should, in principle, also be 
present. Nevertheless, at high tan /3, the genuine SUSY contributions should not be significant and can be 
ignored, while the corrections entering in A;, will almost cancel out when the r'*'r~ decays are considered. 
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Hence, the cross section for 66 — )• $ is rather unstable against scale variation and this 
justifies, a posteriori, the choice of a larger domain of variation with k = 3 in this case, a 
choice that does not appear to be a too extreme one when looking at Fig. O. 
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Figure 15: The scale uncertainty bands of the NLO — >■ <& (left) and the NNLO 66 —J' $ (right) 
cross sections at the IHC at 7 TeV as a function of M$; diiferent values k = 2,3 are used and the 
results are shown when the additional constraint l/n < ^pi/fip < k is imposed or not (marked as 
Kb). In the inserts, the relative deviations (compared to the central cross section values) are shown. 

Let us now turn to the estimation of the uncertainties from the parton densities and 
as- The 90% CL PDF+A^^Pa^ uncertainty, with as(M|) = 0.120 ± 0.002 at NLO for 
gg ^ ^ and as(M|) = 0.1171 ± 0.0014 at NNLO for 66 ^>, is evaluated within the 
MSTW parametrization when including the experimental error on Us- To that, we add in 
quadrature the effect of the theoretical error on a^, estimated by the MSTW collaboration 
to be A*^a, ^ 0.003 at NLO and A*^' ^ 0.002 at NNLO, using the MSTW fixed 
grid with central PDF sets. The 90%CL PDF, PDF+A^^^Pa^ and the PDF+A'^^P+^'^a^ 
uncertainties at the IHC are shown in Fig. ^ as a function of M$. In both the gg ^ ^ and 
66 — )• $ processes, the total uncertainty is below ±10% for M$ < 200 GeV but increases 
at higher masses, in particular in the 66 — )• $ case. 

For completeness, we also display the two cross sections when ones adopts two other 
PDF sets, ABKM and (G) JR, and compare the results with that of MSTW. As can be seen 
in Fig. 17, the deviations from the MSTW values are moderate in the case of 66 — t- <I> for 
the JR scheme: a few percent at M$ = 100 GeV, increasing to ^ 10% at M<j, = 500 GeV. 
The ABKM scheme leads to substantial deviations at high masses, ~ 30% at M$ = 500 
GeV. In the gg fusion process, the prediction with the GJR PDF set at Af$ = 100 GeV 
is ~ 20% lower than in the MSTW and ABKM cases which give comparable results. The 
GJR and ABKM parameterizations cross at ~ 200 GeV and at higher masses, the gg cross 
sections with the ABKM parametrization are ~ 20% lower than for MSTW results while 
the predictions with the GJR set are ~ 10% higher than for MSTW. 

Finally, there is the effect of the uncertainty on the 6-quark mass that had only a 
marginal impact in the SM Higgs case but which will be rather important here. In fact, 
there are three sources of uncertainties which can be attributed to the 6-quark mass. 
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Figure 16: The PDF 90% CL PDF, PDF+A'=''Pa,, and PDF+ A'^'^Pa^ + A^'^a^ uncertainties in the 
MSTW scheme in the — >■ <I> (left) and 66 — > $ (right) cross sections at the IHC at 7 TeV as a 
function of Af$. In the inserts, the relative deviations are shown. 
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Figure 17: The cross sections in the — > $ (left) and 66 $ (right) at the IHC at 7 TeV as a 
function of Af* evaluated when using the ABKM and (G) JR PDF sets. In the inserts, the relative 
deviations from the value in the MSTW scheme are shown. 

The first one is of purely theoretical nature and is due to the choice of the renormaliza- 
tion scheme for the 6-quark mass. In our analysis, we have adopted the MS renormalization 
scheme mainly because of the fact that the calculation of the 66 — )■ <I> process is available 
only in this scheme and we to have chosen to treat the same way both the 66 — t- <I> and 
gg ^ ^ channels. Nevertheless, one could choose another renormalization scheme such as 
the on-shell scheme as it was discussed in the case of the 6-loop contribution in the gg 
fusion process for SM Higgs production. In the MSSM case, one could also adopt the DR 
scheme which appears to be more consistent theoretically. To estimate this scheme depen- 
dence, one can evaluate the difference of the gg ^ ^ cross section in the cases where the 
6-quark mass is defined in the on-shell and in the MS schemes. Using the program HIGLU, 
one obtains a ~ +15% difference. Bearing in mind the fact that the corrections could have 
been negative if we had adopted another scheme (as would have been the case if we have 
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used the DR scheme for instance, although the difference from the MS result would have 
been at the level of a few percent only), one could assign an error A^^^^™^ « ±15% to the 
gg ^ ^ cross section. This is the procedure that we will adopt here. 

Another way to estimate the renormalization scheme dependence of the 6-quark mass, 
would be to look at the differences that one obtains by using rni,{^mh) and mj,{2fni)) as 
inputs in the gg ^ ^ cross section^^ [60]. In this case, one obtains an uncertainty that is 
in fact much larger than the difference between the on-shcll and MS schemes and which 
goes both ways, A^|;<="^<= ^ -20%, +40% for Mq, = 90 GeV for instance. 

In the 66 — > $ process, we cannot perform this exercise as the cross section, using the 
program bbhOnnlo of Ref . [77] , can only be evaluated in the MS scheme with the Yukawa 
couplings evahiatcd at the scale fiR. This is the reason why we extended the domain of 
scale variation in this case to |/xo ^ < 3;Uo- The larger scale uncertainty obtained 

this way could be seen as indirectly taking care of the scheme dependence. 

The second source of uncertainty is of parametric nature and is the same as the one 
affecting the H ^ bb partial decay width of the Higgs boson discussed in section 3. It is 
estimated as previously, i.e. by evaluating the maximal values of the cross sections when 
one includes the error on the input 6-quark mass at the scale nib, rnii{mf,) = 4.19^q[Iq 
GeV, and in the case of the 66 — ^ <1> process where the Yukawa coupling is defined at the 
high scale, the strong coupling constant, ac,.(M^) = 0.1171 ±0.0014 at NNLO, used to run 
the mass mh{fhb) upwards to mbdj^R). In the considered Higgs mass range, one obtains 
an uncertainty of Aj^^^"* Ri —4%, +14% and ^ —3%, +10% in the case of, respectively, 
the gg ^ ^ and 66 ^ $ processes (using the central MSTW PDF set). The difference is 
mainly due to the fact that the bottom quark masses in the two processes are not defined 
at the same scale and, also, in the case of the gg ^ ^ process, additional corrections which 
involve the 6-quark mass, oc log(?fi^/M|), occur at leading order. 

Finally, a third source of uncertainty originates from the choice of the 6-mass value 
in the 6-quark densities. The MSTW collaboration has recently released a set of PDFs 
with different bottom quark masses [55] : it involves six different central PDFs with a range 
of on-shell mi, values between 4.00 GcV and 5.50 GeV in 0.25 GeV steps, in addition to 
the central best-fit with mf, = 4.75 GeV. In order to distinguish between the parametric 
uncertainty in nib and the one due to the correlated PDF-Am^ uncertainty, we have cho- 
sen to calculate the latter uncertainty by taking the minimal and maximal values of the 
production cross sections when using the central value mb = 4.75 GeV and the two closest 
ones upwards and downwards, i.e. nib = 4.5 GeV and = 5 GeV in the MSTW PDF 
set^^. However, we kept in the partonic calculation the central value of mft(mft) = 4.19 GeV 
which, approximately corresponds to the pole mass nib = 4.75 GeV. One obtains a ^ 3-5% 

"^"^Wc should note that the scheme dependence actuaUy appears also as a result of the truncation of the 
perturbative series which, in principle, should already be accounted for by the scale variation. However, the 
scales which enter in the 6-quark mass, that is defined at itself, and in the rest of the gg <b matrix 
element, which is the Higgs mass or the scale jiR, are different. We have checked explictly that the scale 
uncertainties due to the variation of fiR (and /uf) in both schemes are comparable. Adding this scheme 
dependence to the scale variation, as we will do here, is similar in practice to increase the domain of scale 
variation from the central scale while sticking to a given mass renormalization scheme. 

^^Note that when including the Icr errors on WibiWib), while the upper value corresponds to the pole mass 
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uncertainty depending on the considered Higgs mass range. Note that this uncertainty will 
not only affect the cross section in the 66 —t- <^ process in which the 6-densities play the 
major role, but also the one of the gg ^ ^ channel; however, in this case, the change is 
below the percent level and can be safely neglected. 

The effect of these three sources of uncertainties is displayed in Fig. 18 for gg^^ and 
66 — )• ^> as a function of M$. As can be seen, large uncertainties occur, in particular in 
gg^^ where the ~ 15% scheme uncertainty that is absent in bb^^ dominates. 
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Figure 18: The scheme, parametric and PDF uncertainties due to the 6-quark mass in the — > $ 
(left) and 66 — >■ $ (right) cross sections at the IHC at 7 TeV as a function of Af$. In the inserts, 
the relative deviations are shown. 



4.4 The Higgs decay branching fractions 

In the most general case, the decay pattern of the MSSM Higgs particles can be rather 
complicated, in particular for the heavy states. Indeed, besides the standard decays into 
pairs of fermions and gauge bosons, the latter can have mixed decays into gauge and 
Higgs bosons (and the H bosons can decay into hh states) and, if some superparticles 
are light, SUSY decays would also occur. However, for the large values of tan /3 that we 
are interested in here, tan /3 > 10, the couplings of the non-SM like Higgs particles to 
bottom quarks and r leptons are so strongly enhanced and those to top quarks and gauge 
bosons suppressed, that the pattern becomes very simple. To a very good approximation, 
the ^ = A or H{h) bosons will decay almost exclusively into 66 and t~^t~ pairs with 
branching ratios of, respectively, ~ 90% and ~ 10%, with the ti decay channel and the 
decay involving gauge or Higgs bosons suppressed to a level where the branching ratios are 
less than 1%. The CP-even h or H boson, depending on whether we are in the decoupling 
or anti-decoupling regime, will have the same decays as the SM Higgs boson in the mass 
range below MJ^^^ < 135 GeV, decays that have been discussed in section 3 (note that for 
M^^^ ~ 135 GeV, we are in the regime where the decays into 66 and WW are comparable 
and, thus, the uncertainties in the Higgs branching ratios are the largest). 

nib « 5 GeV, the lower value does not correspond to nii, = 4.5 GeV; we will however adopt this smaller 
value to estimate the uncertainty as no other choice is possible within the MSTW set. 
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bb Anil, Aas tot 


TT Ami, Aas tot 


90 


on /in +0.9% +0.2% +0.9% 

»U.W _o^^% -0.2% -0.4% 


Q Rn +3.2% +1.9% +3.8% 

9.DU -s.6%, -1.8% -8.8% 


100 


QH 91 +0.9% +0.2% +1.0% 
^"•^-^ -0.4% -0.2% -0.4% 


Q 7Q +3.3% +1.9% +3.8% 
-8.6% -1.8% -8.8% 


110 


on n/1 +0.9% +0.2% +1.0% 

au.U^ -0.4% -0.2% -0.4% 


Q OR +3.2% +2.0% +3.8% 
-8.6% -1.9% -8.8% 


120 


SQ 88 +1.0% +0.2% +1.0% 
oy.OO _o4% -0.2% -0.4% 


in 19 +3.3% +2.0% +3.8% 
-L'J-J-^ -8.6% -1.9% -8.8% 


130 


an 7/1 +1-0% +0.2% +1.0% 

^'y-'^ -0.4% -0.2% -0.4% 


inOfi +3.2% +2.0% +3.8% 

-LU-ZO _g50/^ _^g% _g7y^ 


140 


an fil +1-0% +0.2% +1.0% 

09. Oi _o.4% -0.2% -0.4% 


in -iQ +3.3% +2.0% +3.8% 

1U.09 _g55;j _^gy^ _g70/„ 


150 


an /L8 +1-0% +0.2% +1.0% 
oy.io _p2% -0.4% 


-ir, ro +3.1% +2.0% +3.7% 

lU.dZ _gg5^^ _2g5,^ _ggO/^ 


160 


aQ '^7 +1-0% +0.2% +1.0% 
_o.4% -0.2% -0.5% 


infi-? +3.2% +2.1% +3.8% 
lU.OO _g5^^ _2o% -s.7% 
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OQ OR +1.0% +0.2% +1.1% 
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in 8/1 +3.2% +2.1% +3.9% 
lU.O^ _8.5% _i.9% _8.7% 
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OQ nc; +1.0% +0.2% +1.1% 
09.Ud -0.4%, -0.2% -0.5% 


-inQc +3.2% +2.0% +3.8% 
lU.9d _g,.5%, -2.0% -8.7% 
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19 t^P. +3.1% +2.2% +3.8% 
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-0.5% -0.3% -0.6% 


197c: +3.1% +2.3% +3.9% 
-L'^-'O -8.3% -2.1% -8.6% 



Table 7: The Higgs decay branching ratio into bb and r+r^ final states (in %) for given Higgs mass 
values (in GeV) with the corresponding individual uncertainties as well as the global uncertainties 
assuming Icr uncertainties on the inputs mb(mfc) and q;s(M|). 

For the evaluation of the theoretical uncertainties in the bb and r+r^ decay branching 
ratios of the states, the analysis of section 3 for the SM Higgs boson can be straightfor- 
wardly extended to the MSSM case. Here, one can ignore the error on the input mass of 
the charm quark (as the decays into cc pairs are strongly suppressed) and consider only the 
uncertainties coming from the two other sources: the inputs mb(m;,) and as{M^). Again, 
the impact of a scale variation in the range ^ fJ- ^ 2M$ is negligibly small. The 

uncertainties on the two branching ratios are displayed in Table ^ for some values of the 
Higgs boson mass, together with the total uncertainties when the individual uncertainties 
resulting from the "Icr" errors on the inputs mh{rnb) and as(M|) are added in quadrature. 

The bb branching ratios of the $ states, BR($ 66) 3m^(M$)/[3mf (M$) + ml], 
slightly decrease with increasing M$ as a result of the higher scale which reduces the 6- 
quark mass mfe(M$), but the total uncertainty is practically constant and amounts to less 
than « 3% as a consequence of the almost complete cancellation of the uncertainties in the 
numerator and denominator. In contrast, there is no such a cancellation in the branching 
fraction for Higgs decays into r+r" pairs, BR($ — )• r+r~) m^/[3m^(M$) + m^], and 
the total uncertainty, that is dominated by the error on the 6-quark mass, reaches the level 
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of 10% for all Higgs masses^*^. If the error on the input 6-quark mass is ignored, the total 
uncertainty in the t^t~ branching ratio will reduce to ~ 3%. 

Finally, let us note that in the MSSM at high tan /3, the total decay widths of the $ 
particles should be taken into account. Indeed, they rise as r($) oc M$ tan^ /3 and thus, 
reach the level of 0(10 GeV) for say, M$ « 200 GeV and tan/3 ^ 50. The total Higgs 
width can thus possibly be larger than the experimental resolution on the r^r^ and bb 
invariant masses when decays into these final states are analyzed. If the total width has to 
be taken into account in the experimental analyses, the uncertainties that affect it should 
also be considered. These uncertainties are in fact simply those affecting the <!'—)• 66 partial 
widths and thus, to a good approximation, the <I> — t- r+r^ branching ratio. The numbers 
given in Table ^ for the uncertainties of BR(<I> — )• r+r^) thus correspond (when multiplied 
by a factor ~ 1.1) to the uncertainties on the total width r(<I>) with a good accuracy. 

4.5 Combined uncertainties 

Let us come to the delicate issue of combining all the uncertainties which come from the 
various sources. Here, we will adopt the procedure A of section 2.3 also used in Ref. [18], 
in which the overall uncertainty on the production cross section is obtained by applying 
the PDF+A^^^as+A^^as uncertainty on the maximal and minimal cross sections from the 
scale variation, adding linearly the scheme/EFT uncertainty. 

As the difference from the SM case is simply the presence of the additional sources 
of uncertainties due to the 6-quark mass, the combination problem reduces in fact to 
answering the question of how one should add these three uncertainties in the production 
cross sections. The uncertainty from the scheme dependence in the gg ^ ^ process should 
be simply, i.e. linearly, added to the scale uncertainty as both emerge as a result of 
the truncation of the perturbative series and are thus of pure theoretical nature. The 
uncertainties of the input 6-quark mass which appears in the parton densities should be 
added to the other PDF+ag uncertainties. As this uncertainty is of experimental origin, 
it should be combined in quadrature with the PDF+A'^^Pas + A^^as uncertainty and 
since it is much smaller than the latter, it will have practically no impact on the total 
PDF+Os error^^. Finally, the parametric uncertainty, which has a special status as it 
appears also in the Higgs branching ratios, can be simply added linearly to the combined 
scale+scheme+PDF uncertainty; we will see below that it will have no impact in practice. 

^^We should note for completeness, that the same results are obtained in the case of the charged Higgs 
particles, where for > rrit + mt ~ 175 GeV and tan/3 ^ 1, one has to consider only the two decay 
modes, -> tb and t~^u with again, BR(//+ tb) ^ 3rnl{M„±) /[3fnl{MH±) + ml] and 

BR(//"'' — ^ Tv) ~ ml /[3rnl{Mjj±) + ml]. The values of the branching ratios, together with the associated 
uncertainties, are thus also those given in Table ^ when AIg> is replaced by Mu± . 

■^^ However, one can also view this uncertainty as being due to the parametrization of the bottom-quark 
PDF which is a theoretical problem, and thus add it linearly to the PDF+Os so that it increases slightly 
to the total PDF uncertainty. We will refrain from doing so, though. Note that one may also consider 
an uncertainty related to the charm-quark mass dependence of the PDFs. This can be estimated again 
as in footnote 9 and we find at most ~ ±2.6% at M$ — 500 GeV in both channels. We then neglect this 
uncertainty added in quadrature to that of the other PDF-related. 
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Applying the procedure above for combining the uncertainties, the results for the 
processes gg ^ ^ and 65 — )• <I> at the IHC with ^/s = 7 TeV are displayed in Fig. as a 
function of M$; the numerical values can also be found in the relevant columns of tables 8 
and 9 where the individual and total uncertainties are given. One finds a total uncertainty 
of +60%, -40% for a{gg $) and « +50%, -35% for a{bb in the low Higgs 

mass range and slightly less at higher masses. As mentioned previously, we expect that 
these numbers approximately hold at least for slightly higher energies, ^/s = 8-10 TeV, 
and probably also at the full fledged LHC with ^/s = 14 TeV. 

Finally, the total uncertainty on the cross section times branching ratio, A(a x BR) 
is obtained by adding the total uncertainties on the production cross sections and the 
uncertainties on the branching fraction in Higgs decays into r+r^ pairs, which is the most 
relevant detection channel. In this addition, at least in the 66 — )• <I> process where one 
defines the $66 coupling at the scale fiR, the uncertainty on the input 6-mass, which is 
common to a and to BR, almost cancels out; only ~ 10% of the error is left out^*^. In the 
case of 55 — 7- <I> where the Yukawa coupling is defined at the scale in contrast to the Higgs 
decay widths, the errors on Os used for the running from mf, to M$ will induce a remaining 
uncertainty. This uncertainty is correlated in a{gg — t- <I>) and BR(<I> — t- t^t~): smaller 
as values lower (j{gg — )■ <5) oc at LO, and also BR(<I' — r+r^) as the resulting ?fife(M$) 
is higher, thus enhancing/reducing the $ — )■ bb/T~^T~ rates. This uncertainty should thus 
be added linearly to the overall scale+scheme+PDF uncertainty of the cross section. 



The impact of this additional uncertainty is shown by the dotted lines of Fig. 19 where 
the uncertainty A(cr x BR) for tau decays is displayed for the two production processes. 
The impact is negligible in the case of 66 fusion and very modest in gg fusion. 
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Figure 19: The total uncertainties due to the scale, PDF and 6-quark mass in the gg — > $ (left) 
and 66 — > $ (right) cross sections at the IHC with 7 TeV as a function of Af$. The dotted lines 
show the uncertainties when those on the branching rations for Higgs decays into r+r^ final states 
are added. In the inserts, the relative deviations are shown. 



^®This is also the case of the A;, SUSY correction which enters both a{gg,bb — >■ $) and BR($ — >■ r^r 
and which thus cancels in the product. This justifies, a posteriori, why we ignored this correction. 



-41 - 







90 


c^o 90+21.7% +7.9% +26.7% +58.9% +53.5% 
0'iZ.ZO_^g 30/^ -8.1% -17.6% -41.7% -41.5% 


95 


414 40+21.1% +7.8% +26.8% +58.3% +52.8% 
°-17.6% -8.0% -17.7% -41.1% -40.9% 


100 


oon 01 +20.4% +7.8% +27.0% +57.7% +52.3% 
"^^"^ -16.9% -7.9% -17.8% -40.5% -40.3% 


105 


orQ 40 + 19.9% +7.8% +27.0% +57.1% +51.7% 
■-^^ '^-16.3% -7.9% -17.9% -40.0% -39.8% 


110 


-iqy (-.1+19.4% +7.7% +27.1% +56.6% +51.3% 
^ -15.9% -7.8% -17.9% -39.6% -39.4% 


115 


1 "i? 'i4+18-9% +7.7% +27.2% +56.2% +51.0% 
-15.4% -7.8% -18.0% -39.2% -39.0% 


120 


19/^ 09+18.5% +7.7% +27.3% +56.1% +50.7% 
"■"^ -15.1% -7.8% -18.0% -39.0%, -38.9% 


125 


^,),, -^j + lS,l^ -7,7^ +27,1^ +o0.1% 
-15.3% -7.7% -IK.IV;. -39.2'/, -39.1% 


130 


00 94+17.7% +7.7% +27.4% +55.1% +50.0% 
-15.1% -7.7% -18.1% -39.0% -38.9% 


135 


fiS 9C + 17.4% +7.7% +27.5% +54.8% +49.6% 
"° "-^-14.7% -7.7% -18.2% -38.7% -38.7% 


140 


c-f? 01+17.1% +7.7% +27.6% +54.5% +49.3% 
dO.Oi_j^4 2% -7.7% -18.2% -38.4% -38.3% 


145 


4f) 7'? + 16-7% +7.6% +27.6% +54.1% +48.9% 
''^-13.8% -7.7% -18.3% -38.1%, -38.2% 


150 


00 QO + 16.4% +7.6% +27.7% +54.0% +48.7% 
'^°-^°-13.5% -7.7% -18.4% -37.8% -37.9% 


155 


09 «S + 16.1% +7.6% +27.8% +53.7% +48.6% 
''^■"°-13.1% -7.7% -18.4% -37.6% -37.5% 


160 


97 c;0+15-9% +7.6% +27.9% +53.5% +48.3% 
' '"-^ -13.0% -7.7% -18.5% -37.5% -37.5% 


165 


90 90+15.7% +7.6% +27.9% +53.3% +48.2% 
'^'^■^^-13.0% -7.7% -18.5% -37.6% -37.5% 


170 


1 q YQ+15.4% +7.6% +28.0% +53.0% +47.9% 
^^•'^-13.0% -7.7% -18.6% -37.7% -37.7% 


175 


1 fi 88+15.2% +7.6% +28.0% +52.8% +47.7% 
■'■"•°°-13.1% -7.7% -18.6% -37.8% -37.7% 


180 


14 4r+15.0% +7.7% +28.1% +52.7% +47.7% 
J-'*-'*J_13.1% -7.7% -18.6% -37.8% -37.8% 


185 


19 49+15.0% +7.7% +28.1% +52.7% +47.5% 
-13.1% -7.7% -18.7% -38.0% -38.2% 


190 


in 71+14-9% +7.6% +28.2% +52.8% +47.6% 
-"--13.2% -7.8% -18.8% -38.1% -38.2% 


195 


q 97+14.9% +7.7% +28.2% +52.7% +47.6% 
-13.2% -7.8% -18.8% -38.2% -38.3% 


200 


^4+14.8% +7.7% +28.0% +52.7% +47.8% 
"•^^-13.2% -7.8% -18.7% -38.3% -38.3% 


225 


4 19 + 14.7% +7.8% +28.3% +52.6% +47.6% 
*--'-^-13.4% -8.0% -19.0% -38.8% -39.0% 


250 


9 94+14.5% +7.9% +28.4% +52.8% +47.8% 
^•^^-13.6% -8.2% -19.2% -39.3% -39.6% 


275 


1 90 + 14.4% +8.0% +28.7% +53.0% +48.3% 
-L-^o_13.7% -8.4% -19.4% -39.8% -39.9% 


300 


r, yf?+14.2% +8.2% +28.9% +53.2% +48.2% 
"•'"-13.8% -8.6% -19.5% -40.3% -40.6% 


325 


r, 47+14.2% +8.4% +29.0% +53.4% +48.5% 
"•^'-13.9% -8.9% -19.6% -40.7% -41.0% 


350 


Q or) + 14.2% +8.7% +29.1% +53.8% +48.9% 
-14.0% -9.2% -19.8% -41.3% -41.7% 


375 


n. iq + 14.1% +9.0% +29.4% +54.2% +49.5% 
^^•-^^-14.1% -9.5% -19.9% -41.8% -41.9% 


400 


1 + 14.1% +9.3% +29.4% +54.5% +49.8% 
"•-^■^-14.2% -9.8% -20.0% -42.1% -42.4% 


425 


r. r|q+14.1% +9.6% +29.7% +55.1% +50.5% 
"•"^-14.3% -10.1% -20.2% -42.6% -42.9% 


450 


r, r)^+14.1% +9.9% +29.8% +55.4% +50.8% 
"•""-14.4% -10.4% -20.3% -43.1% -43.5%, 


500 


r|o + 14.2% +10.6% +30.0% +56.5% +51.9% 
"•""J-14.5% -11.0% -20.5% -44.0% -44.3% 



Table 8: The Higgs production cross sections in the gg ^ ^ channel (for tan/3 = 1) as well as as 
the individual uncertainties (first from scale, then from PDF+ A'^'^p +*'^Q!s at 90%CL and from the 
input mass mi, at Icr) and the total uncertainties for selected values of the Higgs mass. The last 
column displays the total uncertainty when including the combination with the total uncertainty 
on $ — )• T+r~ branching ratio. 
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Mi 


At°t xBR 


90 


487 OQ+26.8% +6.2% +17.8% +49.8% +44.3% 
^°'-°"'-21.1% -8.0% -10.7% -32.8% -32.5% 


95 


400 c;q+25.4% +7.1% +17.7% +47.8% +42.4% 
^"^•''^-19.5% -7.3% -10.3% -31.5% -31.3% 


100 


oAo rn+23.9% +6.4% +17.8% +46.9% +41.4% 
"'^"•"^"-18.4% -7.8% -9.9% -29.9% -29.8% 


105 


904 c7+22.2% +6.3% +18.1% +45.4% +40.0% 
^^^•°'-17.8% -8.3% -9.6% -29.5% -29.3% 


110 


oro rio+21.3% +6.9% +17.2% +43.8% +38.5% 
ZdZ.UO_^g^^^ -7.7% -10.1% -28.9% -28.8% 


115 


91 ^ +6.8% +17.5% +43.0% +37.7% 
^-^^•"^-15.9% -8.0% -9.5% -27.5% -27.3% 


120 


-roe cf; + 18.8% +7.2% +17.6% +41.7% +36.4% 
-^°"-^"-15.4% -7.7% -9.3% -26.9% -27.0% 


125 


Ifil SA+l''-'*'^ +7.7% +17.1% +39.4% +34.1% 
10±.0^_^5 2% -7.6% -9.4% -26.9% -26.8% 


130 


1/in 07+16.4% +7.6% +17.5% +38.8% +33.7% 
■^^^•^'-14.9% -7.8% -9.5% -26.8% -26.7% 


135 


199 0S+-'-^-^'^ +7.3% +17.2% +38.0% +32.8% 
■^^^•^°-14.6% -8.0% -9.2% -26.6% -26.6% 


140 


ins nS+14-1% +''•5% +16-6% +35.6% +30.5% 
■^^°-"°-14.7% -8.0% -9.5% -26.9% -26.8% 


145 


04 fic: + 13.7% +7.7% +17.1% +35.5% +30.2% 
^^■°"^-14.7% -8.2% -9.2% -26.9% -27.0% 


150 


QQ 00 + 13.0% +8.2% +17.0% +35.0% +29.7% 
'5'J-0'J_14.3% _8.0% -9.1% -26.4% -26.5% 


155 


74 1 •3 + 12.0% +8.2% +17.1% +34.0% +28.9% 
'^•-'-'^-14.2% -8.2% -9.0% -26.2% -26.1% 


160 


fjr Q/f + 11.0% +8.3% +16.9% +32.7% +27.6% 
"'J-0'i^l4.1% -8.2% -9.1% -26.4% -26.5% 


165 


ro f!n+10.4% +8.2% +16.9% +32.6% +27.5% 
OO.DU_^3 go/^ -8.5% -8.8% -25.4% -25.5% 


170 


ro 90+9.6% +8.5% +16.6% +31.7% +26.5% 
'J^-^0_13.6% -8.3% -9.2% -25.7% -25.7% 


175 


4fi 7c:+9-l% +8.4% +16.7% +31.3% +26.2% 
^"•'^^-13. 4% -8.6% -8.9% -25.8% -25.8% 


180 


41 07+^-''''^ +8.5% +16.9% +30.4% +25.4% 
'*J--»'_i3.8% -8.7% -8.7% -25.7% -25.7% 


185 


07 CQ+7.5% +9.2% +16.8% +29.8% +24.6% 
<5' -O^-is.s^o _8.2% -8.8% -25.6% -25.8% 


190 


qq 07+7.5% +8.6% +16.5% +29.6% +24.4% 
00.0(_^2.8% -9.2% -8.9% -25.5% -25.6% 


195 


on cc+7.0% +8.9% +16.4% +29.0% +23.9% 
OU.OO_^2.3% -9.1% -9.1% -25.4% -25.5% 


200 


97 fil +6.3% +8.7% +16.4% +28.6% +23.7% 
^'•"-^-12.5% -9.6% -8.8% -25.5% -25.6% 


225 


17 in+^-l^ +9.6% +16.4% +27.6% +22.6% 
-^'■-^"-11.9% -9.7% -8.9% -24.8% -25.0% 


250 


-1 n 07+5.0% +10.6% +16.1% +27.5% +22.5% 
-^"•^'-11.5% -10.0% -8.8% -25.1% -25.4% 


275 


7 97+5.0% +10.9% +16.3% +28.6% +23.9% 
'•^'-10.6% -10.8% -8.6% -24.5% -24.6% 


300 


4 04+^-^''^ +11.6% +16.4% +29.3% +24.3% 
^•^^-10.5% -11.4% -8.4% -24.9% -25.3% 


325 


q ^q+4.8% +12.4% +15.9% +29.7% +24.9% 
^•^•^-9.8% -11.8% -8.6% -24.9% -25.2% 


350 


9 ^q+5.0% +13.1% +16.2% +30.8% +25.9% 
^•^"^-9.4% -12.3% -8.4% -24.6% -25.0% 


375 


1 7r+4.9% +14.1% +16.0% +31.6% +26.9% 
■^•'"^-9.1% -12.5% -8.2% -24.3% -24.5% 


400 


1 90+4.7% +14.9% +15.9% +31.6% +26.9% 
■^•^°_8.9% -13.1% -8.2% -24.8% -25.1% 


425 


n Of.+4-7% +15.5% +16.0% +32.5% +27.8% 
"•^^-8.6% -13.6% -8.3% -25.0% -25.3% 


450 


(1 7-1 +4.7% +16.1% +16.1% +33.4% +28.8% 
U./l_g 2% _14.2% -8.2% -25.2% -25.6% 


500 


n A-, +4.7% +17.8% +15.8% +34.4% +29.8% 
^■^-^-7.9% -15.1% -8.2% -26.1% -26.5% 



Table 9: The same as in Table. 8 but for the bb ^ ^ channel. 
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Thus, the theoretical uncertainties in the bb, gg ^ ^ ^ t~^t~ production and decay 
channels for the MSSM CP-odd and one of the CP-even Higgs particles are extremely 
large. This has important consequences on the MSSM [tan (5, Ma] parameter space that 
can be probed at the IHC. In particular, if no Higgs boson is observed and only exclusion 
limits on the parameter space can be imposed, these uncertainties which can lower the 
production cross sections by a factor of two will have a significant impact. For instance, 
as the cross sections grow as tan^ /3, the values of tan /? which can be ruled out for a given 
A boson mass should be multiplied by a factor fa ^/2 in view of these large uncertainties. 

5. Conclusions 

In the first part of this paper, we have performed a detailed analysis for Standard Model 
Higgs production at the IHC with y/s = 7 TeV in the dominant gluon-gluon fusion chan- 
nel, gg — )• H. We have first updated the production cross sections, including the relevant 
higher order corrections in perturbation theory and then, discussed the theoretical un- 
certainties that affect the predictions: the scale uncertainties that are usually viewed as 
an estimate of the unknown higher order corrections, the uncertainties due to the parton 
distribution functions and the strong coupling constant, and finally, the uncertainties in 
using an effective field theory approach in calculating the higher order corrections beyond 
next-to-leading perturbative order. 

We find that the scale uncertainties, estimated by varying the renormalization and 
factorization scales in the interval ^jUo < I^RiI^f < 2|Lto with ijlq = \Mh being the central 
scale, are moderate, ^ ±12% in the low mass range Mh ~ 120 GeV, and fa +5%, —8% in 
the high mass range, Mh ~ 500 GcV. For the PDF uncertainty, we have considered the 
four NNLO PDF sets that are available, MSTW, ABKM, JR and HERAPDF, and shown 
that there is a sizable spread in the predictions of the cross section. This spread could 
be partly accounted for by taking into account the uncertainties, both experimental and 
theoretical, on the value of ag. The evaluation in the MSTW scheme of the correlated 
PDF+A^'^P+^'^as uncertainties in the gg ^ H cross section leads to a « ±10% uncertainty 
in the entire Higgs mass range that is relevant at the IHC, Mh < 500 GeV. A third source 
of uncertainties is due to the use of the approximation of an infinite mass for the particles 
running in the gg ^ H loop when calculating the amplitude beyond NLO, i.e. for the 
NNLO QCD and the mixed QCD-electroweak corrections; an additional uncertainty arises 
from the renormalization of the bottom quark mass. These uncertainties are quite small 
when taken individually, at most a few percent each, but they add up to a non-negligible 
amount, 3 to 7% depending on the Higgs mass range and should be thus included. 

We have then addressed the issue of combining these three sources of uncertainties. 
Arguing that the PDF uncertainty should be viewed as a pure theoretical one, we have 
proposed three procedures for combining it with those due the scale and the EFT/scheme, 
the simplest one being a linear addition. All three procedures lead to approximately the 
same total uncertainty on the gg ^ H cross section at the IHC, ~ —25%, +30% in the low 
mass range Mh ^ 200 GeV and ~ ±20% at higher Higgs masses. The overall uncertainty 
at the IHC is significantly smaller than the one affecting the cross section at Tevatron 
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energies. This a mere consequence of the smaller scale uncertainty, as the QCD corrections 
are more moderate at IHC energies, and the fact that, at the Tevatron, one is probing the 
high Bjorken-x regime for the gluon densities which is more uncertain. 

We have extended our analysis to other possible energies: ^/s = 8-10 TeV that are 
planed to be explored after the -^/s = 7 TeV run and the LHC design energy of ^/s = 14 
TeV. We provided numerical results for the gg H cross section and in the case of the 
LHC, we also estimated the associated theoretical uncertainties. The later turned out to 
be comparable to those obtained at 7 TeV and these uncertainties are, thus, expected to 
be also the same at the intermediate energies. 

In a second part of the paper, we have addressed the issue of the Higgs decay branching 
ratios and the uncertainties that affect them, namely, the parametric ones due the errors 
in the determination of the bottom and charm quark masses and the QCD coupling a^. 
We obtain uncertainties that are significant in the critical intermediate Higgs mass range 
120 GeV < Mff < 150 GeV where the main Higgs decay channels H ^ bb and H — > WW 
have comparable rates. These uncertainties can reach the level of ^ ±3-10% in the channels 
H WW, ZZ, 77 used to detect the Higgs particle at the IHC; this is also the case for the 
decays H ^ bb and t~^t~ which are also being considered in Tevatron searches. One should 
therefore also take into account these uncertainties which, up to now, have been overlooked 
by the experimental collaborations not only at the LHC but also at the Tevatron. 

Finally, in a third part of the paper, we have extended our analysis of Higgs production 
at the IHC to consider the case of the MSSM neutral Higgs particles. We have investigated 
the production of the CP-odd like particles, ^ = A as well as one of the CP-even H or 
h particles (depending on whether we are in the decoupling or anti-decoupling regimes) 
that arc degenerate in mass with A and have the same couplings, in particular, enhanced 
couplings to bottom quarks for the high values of tan (3 than can be probed at the IHC. 
The two main production processes gg ^ and bb ^ ^, which could have cross sections 
that are order of magnitude larger than in the SM Higgs case, have been considered. 
Numerical results at energies s/s = 7-14 TeV have been given and the associated theoretical 
uncertainties have been evaluated. The latter are due not only to the scale and PDF+a^ 
uncertainties which appear in the SM case, but there are also uncertainties associated to 
the 6-quark which plays a major role in the MSSM. The error on the input 6-quark mass, 
the scheme dependence in the renormalization of the mass of the 6-quark in its contribution 
to the gg ^ ^ amplitude and the effect of nn, on the bottom quark densities in the process 
66 — > $ will induce additional uncertainties. To these, one has to add the parametric 
uncertainties in the branching ratio of the Higgs decay in tau lepton pairs, which is the 
cleanest detection channel at hadron colliders, that is found to be at the level of 10%. 

The overall theoretical uncertainty in these two processes turn out to be extremely 
large at IHC energies, of the order of 50%. This large uncertainty will have a significant 
impact on the MSSM parameter space that can be probed. This is particularly true in the 
case where no MSSM Higgs signal is observed and only exclusion bounds can be derived 
by comparing the experimental data with the predicted Higgs production cross sections. 
As the cross sections in both the gg ^ ^ and 66 — )■ <1> processes increase with tan^ f3, the 
values of tan /3 which could be excluded in the absence of a signal will be smaller by a 
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Figure 20: Same as in Fig. | for = 8 TeV. 

factor ~ \pl if these theoretical uncertainties are taken into account^"^. 
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Al. Addendum: SM cross sections and uncertainties at = 8 TeV 

In this addendum, we summarize the Higgs production cross sections in the SM and in 
the case of gg ^ H we detail the associated uncertainties, at a cm. energy of ^/s = 8 
TeV, following exactly the discussion presented in section 2. In particular we update the 
numbers presented in Table |^ with more precision and a full account of the other main 



processes. The rates for the four production channels are displayed in Fig. ^ and Table 
10. The scale, PDF (using either only the MSTW set or including also the other available 



NNLO PDF sets) and EFT as well as the total uncertainties are shown in Fig. 21 for the 
gg ^ H process and the relevant numbers are given in Table 11. 



This is also the case for the exclusion bounds that have been obtained by the CDF and DO collaborations 
from negative MSSM Higgs searches at the Tevatron [10]. See Ref. [79] for an analysis of this issue. 



Mh 


^NNLO 


NLO 

qq—fu. qq 


„NNLO ^NNLO 
"qq^HW "qq^HZ 




115 


23757.4 


1772.2 


913.9 


508.7 


162.2 


120 


21832.4 


1702.9 


794.1 


446.0 


144.2 


125 


20126.9 


1639.6 


694.8 


392.7 


128.0 


130 


18608.7 


1571.4 


607.7 


346.1 


114.2 


135 


17252.0 


1476.9 


533.5 


306.6 


102.4 


140 


16032.9 


1436.2 


469.6 


272.1 


91.6 


145 


14933.7 


1363.9 


414.8 


241.9 


82.3 


150 


13934.8 


1279.5 


366.7 


215.8 


74.1 


155 


13009.0 


1254.7 


324.3 


192.5 


67.0 


160 


12063.1 


1180.2 


280.8 


169.0 


60.5 


165 


11143.1 


1138.8 


258.7 


156.6 


54.7 


170 


10360.9 


1112.2 


232.1 


141.1 


49.9 


175 


9688.5 


1058.3 


208.2 


127.0 


45.3 


180 


9065.7 


1016.1 


187.7 


114.1 


41.5 


185 


8483.1 


977.8 


171.0 


104.2 


37.9 


190 


7962.6 


944.9 


155.0 


94.7 


34.6 


195 


7509.0 


900.0 


141.0 


86.2 


31.8 


200 


7105.4 


870.3 


128.1 


78.7 


29.4 


210 


6408.1 


807.9 


106.6 


65.7 


24.8 


220 


5823.0 


745.5 


89.2 


55.0 


21.2 


230 


5326.7 


690.4 


75.3 


46.3 


18.3 


240 


4901.8 


638.7 


63.7 


39.2 


16.0 


250 


4536.1 


607.0 


54.3 


33.3 


13.9 


260 


4221.7 


559.7 


46.4 


28.3 


12.2 


270 


3951.5 


525.4 


40.0 


24.2 


10.8 


280 


3721.0 


494.1 


34.5 


20.8 


9.5 


290 


3524.8 


462.8 


29.8 


17.9 


8.5 


300 


3362.2 


433.0 


25.9 


15.5 


7.6 


Table 10: 


Same as in 


Tabic 


with ^/s — 8 


TeV. 
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Figure 21: Same as in Figs. |, |l^, y for = 8 TeV. 
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Mh 


a [pb] 


Scale 


[%] 


PDF+A°7+"^ [%] 


EFT [%] 


Total [%] 


115 


23.76 


+11.6 


-9.6 


+8.9 


-8.5 


±8.7 


+32.3 


-25.0 


120 


21.83 


+12.8 


-9.5 


+8.9 


-8.5 


±8.7 


+35.0 


-24.9 


125 


20.13 


+11.9 


-9.5 


+8.9 


-8.5 


±8.7 


+35.3 


-24.8 


130 


18.61 


+10.3 


-9.4 


+8.9 


-8.5 


±8.7 


+32.7 


-24.7 


135 


17.25 


+9.4 


-9.4 


+8.9 


-8.4 


±8.6 


+30.0 


-24.6 


140 


16.03 


+8.9 


-9.3 


+8.9 


-8.4 


±8.6 


+28.5 


-24.5 


145 


14.93 


+8.5 


-9.3 


+8.9 


-8.4 


±8.5 


+27.6 


-24.4 


150 


13.93 


+8.3 


-9.2 


+8.8 


-8.4 


±8.4 


+27.0 


-24.3 


155 


13.01 


+8.1 


-9.2 


+8.8 


-8.5 


±8.3 


+26.4 


-24.1 


160 


12.06 


+8.0 


-9.1 


+8.8 


-8.5 


±7.6 


+25.5 


-23.4 


165 


11.14 


+7.8 


-9.1 


+8.8 


-8.5 


±6.6 


+24.3 


-22.5 


170 


10.36 


+7.7 


-9.1 


+8.8 


-8.5 


±5.9 


+23.3 


-21.7 


175 


9.69 


+7.6 


-9.0 


+8.9 


-8.5 


±5.4 


+22.7 


-21.2 


180 


9.07 


+7.5 


-9.0 


+8.9 


-8.5 


±4.9 


+22.0 


-20.6 


185 


8.48 


+7.5 


-9.0 


+8.9 


-8.5 


±4.7 


+21.9 


-20.5 


190 


7.96 


+7.5 


-9.0 


+8.9 


-8.5 


±5.1 


+22.2 


-20.9 


195 


7.51 


+7.4 


-9.0 


+8.9 


-8.6 


±5.3 


+22.4 


-21.1 


200 


7.11 


+7.4 


-8.9 


+8.9 


-8.6 


±5.4 


+22.4 


-21.2 


210 


6.41 


+7.3 


-8.9 


+8.9 


-8.7 


±5.4 


+22.4 


-21.2 


220 


5.82 


+7.2 


-8.9 


+8.9 


-8.7 


±5.3 


+22.1 


-21.1 


230 


5.33 


+7.0 


-8.9 


+9.0 


-8.8 


±5.1 


+21.8 


-20.9 


240 


4.90 


+6.9 


-8.8 


+9.0 


-8.8 


±4.9 


+21.4 


-20.7 


250 


4.54 


+6.7 


-8.8 


+9.1 


-8.9 


±4.6 


+21.0 


-20.5 


260 


4.22 


+6.7 


-8.8 


+9.1 


-8.9 


±4.4 


+20.7 


-20.4 


270 


3.95 


+6.6 


-8.8 


+9.2 


-8.9 


±4.1 


+20.5 


-20.2 


280 


3.72 


+6.6 


-8.8 


+9.3 


-9.0 


±3.9 


+20.3 


-20.0 


290 


3.52 


+6.5 


-8.7 


+9.3 


-9.1 


±3.7 


+20.1 


-19.8 


300 


3.36 


+6.5 


-8.7 


+9.3 


-9.2 


±3.4 


+20.0 


-19.7 


320 


3.14 


+6.5 


-8.7 


+9.4 


-9.4 


±2.9 


+19.9 


-19.7 


340 


3.10 


+6.6 


-8.9 


+9.6 


-9.5 


±3.3 


+20.5 


-20.4 


360 


3.33 


+6.7 


-8.9 


+9.8 


-9.6 


±4.1 


+21.4 


-21.4 


380 


3.19 


+6.1 


-8.6 


+10.0 


-9.6 


±3.2 


+20.2 


-20.4 


400 


2.86 


+5.7 


-8.3 


+10.2 


-9.8 


±3.8 


+20.5 


-20.8 


450 


1.94 


+4.9 


-7.9 


+10.6 


-10.1 


±4.6 


+21.1 


-21.5 


500 


1.24 


+4.0 


-7.6 


+11.2 


-10.5 


±5.1 


+21.1 


-22.0 



Table 11: Same as in Table ^ with ^/s = 8 TeV, including only procedure A described in the text 
to calculate the total uncertainty. 



-49- 



References 

[1] P. Higgs, Phys. Lett. 12 (1964) 132; F. Englert and R. Brout, Phys. Rev. Lett. 13 (1964) 
321; G. Guralnik, C. Hagen and T. Kibble, Phys. Rev. Lett. 13 (1964) 585; P. Higgs. Phys. 
Rev. 145 (1966) 1156. 

[2] For a review, sec: A. Djouadi, Phys. Rcpt. 457 (2008) 1 [hcp-ph/0503172]. 

[3] For a rewiew see, M. Drees, R.M. Godbole and P. Roy, Theory and Phenomenology of 
Sparticles, World Scientific, Spring 2004. 

[4] J. Gunion, H. Haber, G. Kane and S. Dawson, The Higgs Hunter's Guide, Addison-Wesley, 
Reading 1990. 

[5] A. Djouadi, Phys. Rep. 459 (2008) 1 [hep-ph/0503173]. 
[6] M. Spira, Fortschr. Phys. 46 (1998) 203. 

[7] S. Heinemeyer, W. HolUk and G. Weiglein, Phys. Rept. 425 (2006) 265; S. Heinemeyer, Int. 
J. Mod. Phys A21 (2006) 2659, B. Allanach et al., JHEP 0409 (2004) 044. 

[8] The LEP GoUaborations, Phys. Lett. B565 (2003) 61. 

[9] The CDF and DO coll., "Combination of Tevatron Searches for the Standard Model Higgs 
Boson in the W+W- Decay Mode", Phys. Rev. Lett. 104 (2010) 061802; The CDF and DO 
coll., "Combined CDF and DO upper limits on Standard Model Higgs boson production 
with up to 6.7 fb-i of Data", arXiv: 1007.4587 [hep-ex]. 

[10] The TevNPHWG for the GDF/DO Collaborations, arXiv:1003.3363 [hep-ex]. 

[11] See for instance, R. Heuer, talk given at IWLC2010, CERN Geneva, October 2010. 

[12] ATLAS Collaboration, "ATLAS sensitivity prospects for Higgs boson production at the 
LHC running at 7,8 or 9 TeV", ATLAS-PHYS-PUB-2010-015. 

[13] CMS Collaboration, "The CMS Physics reach for searches at 7 TeV", CMS Note 2010/008. 

[14] LHC Higgs cross section Working Group, "Handbook of LHC Higgs Cross Sections" , in 
preparation; C. Mariotti, S. Dittmaier, G. Passarino and R. Tanaka (eds). 

[15] H. Georgi, S. Glashow, M. Machacek and D. Nanopoulos, Phys. Rev. Lett. 40 (1978) 692. 

[16] M. Dittmar and H. Dreiner, Phys. Rev. D55 (1997) 167. 

[17] S.L. Glashow, D.V. Nanopoulos and A. Yildiz, Phys. Rev. D18 (1978) 1724. 

[18] J. Baglio and A. Djouadi, JHEP 1010 (2010) 064. 

[19] A. Djouadi, M. Spira and P. Zerwas, Phys. Lett. B264 (1991) 440; S. Dawson, Nucl. Phys. 
B359 (1991) 283. 

[20] M. Spira, A. Djouadi, D. Graudenz and P.M. Zerwas, Nucl. Phys. B453 (1995) 17. 

[21] R.V. Harlander and W. Kilgorc, Phys. Rev. Lett. 88 (2002) 201801; C. Anastasiou and K. 
Melnikov, Nucl. Phys. B646 (2002) 220; V. Ravindran, J. Smith and W.L. Van Neerven, 
Nucl. Phys. B665 (2003) 325. 

[22] S. Catani, D. de Florian, M. Grazzini and P. Nason, JHEP 0307 (2003) 028. 

[23] A. Djouadi, M. Spira and P.M. Zerwas, Z. Phys. C70 (1996) 427. 

[24] D. Dicus and S. Willenbrock, Phys. Rev. D39 (1989) 751. 



-50- 



[25] J. Campbell, R. K. Ellis, F. Maltoiii and S. Willcnbrock, Phys. Rev. D67 (2003) 095002; F. 
Maltoni, Z. Sullivan and S. Willenbrock Phys. Rev. D67 (2003) 093005. 

[26] R. Harlandcr and W. Kilgorc, Phys. Rev. D68 (2003) 013001. 

[27] C. Anastasiou, R. Boughezal and F. Pietriello, JHEP 0904 (2009) 003. 

[28] D. de Florian and G. Grazzini, Phys. Lett. B674 (2009) 291. 

[29] S. Moch and A. Vogt, Phys. Lett. B631 (2005) 48; V. Ravindran, Nucl. Phys. B752 (2006) 
173; V. Ahrens, T. Becher, M. Neubert and L.L. Yang, Eur. Phys. J. C62 (2009) 333. 

[30] See for instance, G. Corceha and L. Magnea, Phys. Rev. D72 (2005) 074017. 

[31] See for instance, V. Ravindran, J. Smith and W.L. van Neerven, Nucl. Phys. B634 (2002) 
247; C. Anastasiou, K. Melnikov and F. Petriello, Nucl. Phys. B724 (2005) 197; S. Catani 
and M. Grazzini, Phys. Rev. Lett. 98 (2007) 222002; C. Anastasiou, S. Bucherer and Z. 
Kunszt, JHEP 0910 (2009) 068. 

[32] C. Anastasiou et al., JHEP 0908 (2009) 099. 

[33] C. Anastasiou, "Higgs production via gluon fusion", talk given at Higgs Hunting Workshop, 
Orsay-Prance, 29-31 July 2010 

[34] S. Actis, G. Passarino, C. Sturm and S. Uccirati, Nucl. Phys. B811 (2009) 182. 

[35] For the electroweak corrections, see also: A. Djouadi and P. Gambino, Phys. Rev. Lett. 73 
(1994) 2528; U. Aglietti, R. Bonciani, G. Dcgrassi and A. Vicini, Phys. Lett. B595 (2004) 
432; G. Degrassi and F. Maltoni, Phys. Lett. B600 (2004) 255; S. Actis et al, Phys. Lett. 
B670 (2008) 12. 

[36] E. Berger, C. Qing-Hong, C. Jackson and G. Shaughnessy, Phys. Rev. D82 (2010) 053003; 
F. Demartin, S. Forte, E. Mariani, J. Rojo and A. Vicini, Phys. Rev. D82 (2010) 014002; V. 

Ahrens, T. Becher, M. Neubert and L.L. Yang, arXiv:1008.3162 [hep-ph]; S. Alekhin, J. 
Blumlein, P. Jimenez-Delgado, S. Moch and E. Reya, arXiv:1011.6259 [hep-ph]. 

[37] The Fortran codes can be found in Michael Spira's web page, 

http : //people .web .psi . ch/~ mspira/. For HIGLU, see hep-ph/9510347. 

[38] A.D. Martin, W. Strirling, R. Thorne and G. Watt, Eur. Phys. J. C63 (2009) 189. 

[39] G. Altarelh, R.K. Ellis and G. Martinelh, Nuc. Phys. B157 (1979) 461; J. Kubar-Andre and 
F. Paige, Phys. Rev. D19 (1979) 221; T. Han and S. Willenbrock, Phys. Lett. B273 (1991) 
167; J. Ohnemus and W. J. Stirling, Phys. Rev. D47 (1993) 2722; A. Djouadi and M. Spira, 
Phys. Rev. D62 (2000) 014004. 

[40] R. Hamberg, W.L. van Necrvcn and T. Matsuura, Nucl. Phys. B359 (1991) 343 and (E) 

ibid. B644 (2002) 403; O. Brein, A. Djouadi and R. Harlandcr, Phys. Lett. B579 (2004) 149. 

[41] M. L. Ciccolini, S. Dittmaier and M. Kramer, Phys. Rev. D68 (2003) 073003. 

[42] T. Han, G. Valencia and S. Willenbrock, Phys. Rev. Lett. 69 (1992) 3274; T. Figy, C. Oleari 
and D. Zeppcnfcld. Phys. Rev. D68 (2003) 073005. 

[43] P. Bolzoni, F. Maltoni, S. Moch and M. Zaro, Phys. Rev. Lett. 105 (2010) 011801 and 
arXiv:1006.2323 [hep-ph]. 

[44] M. Ciccolini, A. Denner and S. Dittmaier, Phys. Rev. D77 (2008) 013002. 



- 51 - 



[45] W. Bccnakkcr ct al., Phys. Rev. Lett. 87 (2001) 201805; Nucl. Phys. B653 (2003) 151; S. 
Dawson et al., Phys. Rev. Lett. 87 (2001) 201804 and Phys. Rev. D67 (2003) 071503. 

[46] We are grateful! to Michael Spira for a discussion on this issue. 

[47] R. Harlander and K. Ozeren, arXiv:0909.3420; A. Pak, M. Rogal and M. Steinhauser, 
arXiv:0911.4662; R. Harlander, H. Mantler, S. Marzani, K. Ozeren, arXiv:0912.2104; S. 
Marzani et al. Nucl. Phys. B800 (2008) 127. 

[48] A.D. Martin, W. Strirling, R. Thorne and G. Watt, Eur. Phys. J. C64 (2009) 653. 

[49] P.M. Nadolsky et al. (CTEQ coll.), Phys. Rev. D78 (2008) 013004. 

[50] R.D. Ball et al., Nucl. Phys. B823 (2009) 195. 

[51] P. Jimenez-Delgado and E. Reya, Phys. Rev. D80 (2009) 114011. 

[52] S. Alekhin, J. Blumlein, S. Klein and S. Moch, arXiv;0908.2766. 

[53] The NNLO PDF sets can be found at; www.desy .d.e/hlzeus/combined_results. 

[54] K. Nakamura et al.. Particle Data Group, J. Phys. G37 (2010) 075021. 

[55] A.D. Martin, W.J. Stirling, R.S. Thorne and G. Watt, arXiv:1007.2624. 

[56] PDF4LHC recommandation can be found at: 

http : / / www . hep . ucl . ac. uk/pdf 41hc/PDF4LHCrecom . pdf . 

[57] M. Cacciari, S. Frixionc, M.L. Mangano, P. Nason and G. Ridolfi, JHEP 0809 (2008) 127. 

[58] N. Gray, D.J. Broadhurst, W. Grafe and K. Schilcher, Z. Phys. C48 (1990) 673; K. 
Chetyrkin and M. Steinhauser, Nucl. Phys. B573 (2000) 617; K. Melnikov and T. van 
Ritbergen, Phys. Lett. B482 (2000) 99. 

[59] S. Gorishny, A. Kataev, S. Larin and L. Surguladze, Phys. Rev. D43 (1991) 1633; K. 
Chetyrkin, Phys. Lett. B404 (1997) 161; J. Vermaseren, S. Larin and T. van Ritbergen, 
Phys. Lett. B405 (1997) 327. 

[60] We thank Michael Spira for a discussion on this point. 

[61] K. Chctyrkhin, B. Kniehl and M. Steinhauser, Phys. Rev. Lett. 79 (1997) 353; P. Baikalov 
and K. Chetyrkhin, Phys. Rev. Lett. 97 (2006) 061803. 

[62] A. Bredenstein ct al., Phys. Rev. D74 (2006) 013004 and JHEP 070 (2007) 080. 

[63] A. Ghinculov, Phys. Lett. B337 (1994) 137; L. Durand. K. Riesselnian and B. Kniehl, Phys. 
Rev. Lett. 72 (1994) 2534; V. Borodulin and G. Jikia, Phys. Lett. B391 (1997) 434. 

[64] A. Djouadi, M. Spira, J. van de Bij and P.M. Zerwas, Phys. Lett. B257 (1991) 187; A. 
Djouadi, M. Spira and P.M. Zerwas, Phys. Lett. B311 (1993) 255. 

[65] K. Chctyrkhin ct al., Phys. Rev. D80 (2009) 074010. 

[66] A. Djouadi, J. Kalinowski and M. Spira, Coniput. Phys. Commun. 108 (1998) 56. 

[67] A. Djouadi, N. Kauer and M. Kramer, in progress. 

[68] M. Carena, S. Heinemeyer, C. Wagner and G. Weiglein, Eur. J. Phys. C26 (2003) 601. 

[69] E. Boos et al., Phys. Rev. D66 (2002) 055004; E. Boos, A. Djouadi and A. Nikitenko, Phys. 
Lett. B578 (2004) 384. 



-52- 



[70] D.M. Pierce, J.A. Bagger, K.T. Matchev and R.J. Zhang, Nucl. Phys. B491 (1997) 3. 

[71] See e.g., M. Carena, D. Garcia, U. Nierste and C.E. Wagner, Nucl. Phys. B577 (2000) 88. 

[72] D. Noth and M. Spira, Phys. Rev. Lett. 101 (2008) 181801. 

[73] S. Heinemeyer, W. Hollik and G. Weiglein, Comp. Phys. Commun. 124 (2000) 76. 

[74] A. Djouadi, J.L. Kneur and G. Moultaka, Comput. Phys. Commun. 176 (2007) 426. 

[75] S. Dawson, A. Djouadi and M. Spira, Phys. Rev. Lett. 77 (1996) 16; R. Harlandcr and M. 
Steinhauser, JHEP 0409 (2004) 066, ibid. Phys. Rev. D68 (2003) 111701; C. Anastasiou et 
al., JHEP 0701 (2007) 082; U. Aglietti et al, JHEP 0701 (2007) 021; R. Bonciani, 
G. Degrassi and A. Vicini, JHEP 0711 (2007) 095; M. Muhlleitner and M. Spira, Nucl. Phys. 
B790 (2008) 1; G. Degrassi and P. Slavich, Nucl. Phys. B805 (2008) 267; S. Dawson and C. 
Jackson, Phys. Rev. D77 (2008) 015019; C. Anastasiou, S. Beerh and A. Daleo, Phys. Rev. 
Lett. 100 (2008) 241806; M. Muhlleitner, H. Rzehak and M. Spira, arXiv:101.3214 [hep-phj. 

[76] S. Dittmaier, M. Kramer and M. Spira, Phys. Rev. D70 (2004) 074010; S. Dawson, C. 
Jackson, L. Reina and D. Wackeroth, Phys. Rev. D69 (2004) 074027. 

[77] We thank Robert Harlander for providing us with his code. 

[78] J. Campbell et al., hep-ph/0406152. 

[79] J. Baglio and A. Djouadi, Revisiting the constraints on the Supersymmetric Higgs sector at 
the Tevatron, arXiv:1012.2748. 



-53- 



